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An effective quark Lagrangian is derived from first principles througli bilocal gluon field corre- 
lators. It is used to write down equations for baryons, containing both perturbative and nonper- 
turbative fields. As a result one obtains magnetic moments of octet and decuplet baryons without 
introduction of constituent quark masses and using only string tension as an input. Magnetic mo- 
ments come out on average in reasonable agreement with experiment, except for nucleons and E~. 
The predictions for the proton and neutron are shown to be in close agreement with the empirical 
. values once we choose the string tension such to yield the proper nucleon mass. Pionic corrections 

' to the nucleon magnetic moments have been estimated. In particular, the total result of the two- 

^Sj I body current contributions are found to be small. Inclusion of the anomalous magnetic moment 

^ ■ contributions from pion and kaon loops leads to an improvement of the predictions. 

o 

^ ■ I. INTRODUCTION 

The QCD dynamics of qq and 3g systems is governed by two basic phenomena: confinement and chiral symmetry 
I , breaking (CSB), which should be treated in a fully relativistically covariant way. Confinement is usually introduced 
' for static quarks via the area law of the Wilson loop [1] or equivalently through the field correlators in the Field 

Correlator Method (PCM) [2,3]. 
' For spinless quarks, or neglecting spin-dependent mass corrections, one can envisage a self-consistent method which 
CO ' treats confinement as the area law also for light quarks in a relativistically covariant way. Such method was introduced 
] . originally in [4] for mesons, in [5] for baryons, and in [6] for heavy-light mesons, and later on in [7] the method was 
' I ' generalized taking into account the dynamical degrees of freedom of the QCD string, which naturally appears due to 
^ . the area law. 

As a result Regge trajectories have been found in [7] with the correct string slope {2n(j)^^ . It was realized later on, 
that the method used in [4]- [7] can be more generally developed in the framework of the so-called einbein formalism, 
see [8]- [10]. Spin corrections have been considered in [11] for heavy mesons and in [6] for heavy-light ones. In the 
O il general case of light quarks spin-dependent correlations have been introduced in [12], and for gluons in [13]. For a 
1^ I general review with explicit formulae see [14]. Baryon Regge trajectories have been found in [5]. In all cases the basic 
formalism is the FCM and the Feynman-Schwinger (or world-line) path integral representation [3,15,16] which is well 
suited for relativistic quarks when spin is considered as a perturbation. 

The main difficulty which was always present in the method, was the perturbative treatment of spin degrees of 
freedom (which is incorrect, e.g., for the pion) and absence of spontaneous CSB effects in general [17]. Recently a new 
■ type of formalism was suggested to treat simultaneously confinement and CSB and a nonlinear equation was derived 
for a light quark in the field of heavy antiquark [18]. This equation derived directly from QCD Lagrangian was found 
to produce linear confinement and CSB for the light quark and the explicit form of the effective quark mass operator 
M{x,y) was defined obeying both these properties. 

The eigenvalues and eigenfunctions of the nonlocal and nonlinear equations have been determined and a nonzero 
condensate was computed in [19], confirming that CSB is really present in the equations. In an additional study [20] it 
was demonstrated that magnetic field correlators do not contribute to the large distance confinement, however strongly 
modify the confinement for lowest levels and heavy-light masses corrected in this way are favourably compared in [20] 
to the experiment and results of other calculations. 

Moreover, it was shown in [21] that lattice data strongly support the dominance of the Gaussian (bilocal) correlator, 
estimating the correction due to higher correlators to 1-2%. Since the method of [18] is quite general and allows to 
treat also multiquark systems, it can be applied to the qq and 3q systems, to find dynamical equations for them, 
which contain confinement and CSB [22]. To make these equations tractable, one systematically exploits the large Nc 
limit, and mostly confine ourselves to the simplest field correlators - the so-called Gaussian approximation; it was in 
particular shown in [18] that the sum over all correlators does not change the qualitative results. However, the kernel 
of equations becomes much more complicated. 

In the present paper we study the baryon magnetic moments based on the derived effective Lagrangian without 
constituent quark masses. The paper is organized as follows. In Section 2 the general effective quark Lagrangian 
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from the standard QCD Lagrangian is obtained by integrating out gluonic degrees of freedom, and the nonhnear 
equation for the single quark propagator S (attached to the string in a gauge invariant way) is derived, following the 
procedure in [22]. Section 3 is devoted to the baryon Green's function, which can be expressed in the lowest order of 
our approximation scheme (neglecting gluon and pion exchanges) in terms of 3 independent quark Green's function, 
resulting in a Hamiltonian as a sum of three quark terms. In section 4 the next order approximation is written down 
when perturbative gluon exchanges are taken into account, including the nonperturbative interaction between quarks 
violating the factorized form of the zeroth order approximation. The next section is devoted to the calculation of 
magnetic moments of baryons both in octet and decuplet representations of SU(3) flavour group. In section 6 we 
discuss the corrections to magnetic moments due to pion exchange contributions. 

II. EFFECTIVE QUARK LAGRANGIAN 

As was discussed in the previous section, one can obtain an effective quark Lagrangian by averaging over background 

gluonic fields. We shall repeat this procedure following [18] now paying special attention to the dependence on the 
contour in the definition of contour gauge, and introducing the operation of averaging over contour manifold. The 
QCD partition function for quarks and gluons can be written as 

Z = j DAD^Dilj+exp[Lo + Li + Li„t], (1) 

where we are using Euclidean metric and define 

Lo = -\JdMF;j', (2) 

Li = -ij ^i}+{x){d + ms) ^i){x)(fx, (3) 



X. (4) 

Here and in what follows ^ ijjaa denotes quark operator with flavour /, color a and bispinor index a. 

To express A^{x) through F^i, one can use the generalized Fock-Schwinger gauge [23] with the contour C{x) from 
the point a; to a;o, which can also lie at infinity. 

Now one can integrate out the gluonic field A^(a;), and introduce an arbitrary integration over the set of contours 
C{x) with the weight Dk,{C), since Z is gauge invariant it does not depend on the contour C{x). One obtains 

Z = j DK{C)D^D^+exp{Li + Leff}, (6) 

where the effective quark Lagrangian L^s is defined as 

expLeff = (exp j ^ip+A f'ij)d^x)A- (7) 

Using the cluster expansion. Loff can be written as an infinite sum containing averages {{A)^)a- At this point one 
can exploit the Gaussian approximation, neglecting all correlators {{A)^) of degree higher than k = 2. Numerical 
accuracy of this approximation was discussed and tested in [21]. One expects that for static quarks corrections to 
Gaussian approximation amount to less than 2-3%. 
The resulting effective Lagrangian is quartic in if), 

= 2^ / d'xd^f^tai^) ^Mx) O^liv) 'My)Ja0;^6{x, v) + ©(V'^), (8) 
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and J^v is expressed as 

J,.{x,y)=g'J^ ^du. ^du/^^{F,^{u)F,,^,{v)). (10) 

Leffi Eq. (8), is written in the contour gauge [23]. It can be identically rewritten in the gauge-invariant form if one 
substitutes parallel transporters $(a;,a:o), $(y,a;o) (identically equal to unity in this gauge) into Eqs. (8) and (10), 
multiplying each i/)(a;) and V'(2/) respectively and replacing F{u) in Eq. (10) by <^{x,u)F{u)'^{u,Xq) and similarly for 
F{v). 

After that Leff becomes gauge-invariant, but in general contour-dependent, if one keeps only the quartic term (8), 
and neglects all higher terms. A similar problem occurs in the cluster expansion of Wilson loop, when one keeps only 
lowest correlators, leading to the (erroneous) surface dependence of the result. The situation here is the same as with 
a sum of QCD perturbation series, which depends on the normalization mass ji for any finite number of terms in the 
series. This unphysical dependence is usually treated by fixing at some physically reasonable value /^o (of the order 
of the inverse size of the system) . 

The integration over contours Dk{C) in (6) resolves this difficulty in a similar way. Namely, the partition function Z 
formally does not depend on contours (since it is integrated over a set of contours) but depends on the weight Dk{C). 
We choose this weight in such a way, that the contours would generate the string of minimal length between q and q. 
Thus the physical choice of the contour corresponds to the minimization of the meson (baryon) mass over the class of 
strings, in the same way as the choice oi fi — fio corresponds to the minimization of the dropped higher pcrturbativc 
terms. As a practical outcome, we shall keep the integral Dk{C) till the end and finally use it to minimize the string 
between the quarks. 

Till this point we have made only one approximation -neglected all field correlators except the Gaussian one. Recent 
lattice calculations (see Refs. [24,25]) estimate the accuracy of this approximation at the level of few perccnts. Now 
one must use another approximation, i.e. assume a large Nc expansion and keep the lowest term. As was shown in 
[18] this enables one to replace in (8) the colorless product ^ipb{.x) ^tp^{y) = tr{ ^ tp{x)^{x, xo)^{xo,y) ^tp~^{y)) by 
the quark Green's function 



Mx) (2/) - 6fgN,Si3^{x,y). (11) 



ig^^ assumes the form 



L^^ =-ij d^xd^y fi^tai^) f'M^six, y) ^ Vai(y), (12) 
where the quark mass operator is 

^M^s{x,y) = -J^v{x,y){'y^ ^ S{x,y)j^)a5- (13) 
Prom (12) it is evident that ■''5' satisfies 

{-id^-inif) ^'S{x,y)-i [ ^' M{x,z)d*z ^ S{z,y) = S^^\x - y). (14) 



Eqs. (13)-(14) have been first derived in [18]. From (6) and (12) one should expect that at large Nc the qq and 
3q dynamics is expressed through the quark mass operator (13), which should contain both confinement and CSB. 
Indeed, the analysis performed in Refs. [18]- [20] reveals that confinement is present in the long-distance form of 
M{x,y), when both distances |x|, |y| of light quark from heavy antiquark (placed at x = 0) are large. 

We shall now make several simplifying assumptions, to clarify the structure of M(x,y). First of all we take the 
class of contours C going from any point x = (2:4, x) to the point (0:4,0) and then to (—00) along the X4 axis. For 
this class the corresponding gauge was studied in [26]. Secondly, we take the dominant part of J^^ in (13), namely 
J44, which is proportional to the correlator of color-electric fields. This yields a linear confining interaction, while 
the other components Jik, Ji4, J4i,i = 1,2,3 have been neglected, containing magnetic fields and yielding momentum 
dependent corrections. (It is easy to take into account these contributions in a more detailed analysis). 

The correlator {FF) in (10) can be expressed through the scalar correlator D(x), defined as [2], 

-^{F„p{u)^{u,v)F^sivMv,u)) = D{u - v){5c,^50s - J^^^/Jt) + 0{D^), (15) 
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where the correlator Di, not contributing to confinement, is neglected. As a result one has for M [19,20] 

fMc^^{x,y)= ■fMW/+ ■fM«CTi+ /m(4)74+ •^MW7i. (16) 

Here we have defined 

The dominant part of M, ^M^^^ is linearly growing at large |x|, |y| and in the most simple case of Gaussian form of 
D{x), can be written as 

1 (^4-H4) -v- _l_ V 

fM('\x,y) = ^r^e a\^\~5('\^ - y) (18) 

where Tg is the gluon correlation length, and 5 is a smeared ^-function, which can be represented as [19,20] 

5(3)(x-y)«earp(-t5^)(_I=)3, ^ 2T,. (19) 

Here again T„ is the gluon correlation length, which enters D{u) as D{u) ~ L)(0)exp(— ^^j-). We are now in the 

a 

position to derive the qq, 3g Green's function, which will be done in the next section. 



III. EQUATIONS FOR THE BARYON GREEN'S FUNCTION 

Equations for the 3q system can be written in the same way as for the qq system. We again shall assume the large 
Nc limit in the sense, that 1/Nc corrections from qq pairs to the quark Green's function and the effective mass can 
be neglected. We now write down the explicit expressions for Nc = 3. 
The initial and final field operators are 

^ra{x,y,z) = e„5cr"'^^^aa(2;,C(,x))V^5,fl(?y.C(?y))V^e,,(z,C(z)) (20) 

with the notations: a, b, c, are color indices, a, j3, 7 arc Lorentz bispinor indices and transported quark operators are 

i^accix, C{x)) = {^c{x, X)ll>c,{x))a. (21) 

The contour C{x) in can be arbitrary , but it is convenient to choose it in the same class of contours that is used 
in Dk{C) and in the generalized Fock-Schwinger gauge [23]. T°'l^^ is the Lorentz spinor tensor securing proper baryon 
quantum numbers. One can also choose other operators, but it does not influence the resulting equations. In Eq. (20) 

we have omitted flavour indices in T and ipi^, C), to be easily restored in final expressions. 

Using now the effective Lagrangian (12) valid at large we obtain for the 3g Green's function. 

G^^^\x,y,z\x',y',z') = 

^ J Dn{C)D^I^Di,+ ^f,^{x', y', z')-^Ux, y, z)exp(Li + L^e). (22) 

Integrating out quark degrees of freedom and neglecting the determinant at large Nc, one has 

G(39)= J DK{C){eT){e'T'){S{x,x')S{y,y')S{z,z')+pej:m}, (23) 

where for simplicity color and bispinor indices are suppressed together with parallel transporters in initial and final 
states. 

One can also define unprojected (without F, F') 3q Green's function Gim"^ with 3 initial and 3 final bispinor indices 
instead of projected by F, F' quantum numbers of baryon. Assuming that minimization over contours Dk{C) reduces 
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to the single choice of the contours (the single string junction trajectory minimizing the mass of baryon), one can 
write equation for G^n''* : 



{—idx — im\ — iMi){—idy — im2 — iM2){—idz — imz — iM3)G^ 
X - x')S^^\y - y')5'^^\z - z') 



3q) 
un 



(24) 



with e.g. MiG = J M{x,u)G{u,x')d'^u. One can simphiy the form (23) for G^^^) taking into account that M{x,x') 
actually does not depend on Hence the interaction kernel of G^ '''^ does not depend on relative energies, as 

in [27]. Similarly to [27,28] one can introduce Fourier transform of G*^'^*-' in time components and take into account 
energy conservation E = Ei + E2 + E^. One obtains 



where we have used the notation 



G^^''\E,E2,E:i):^ j DK{C){eT){e'T') 



[E — E2 — Es — Hi)(E2 — H2){Ez — H3) ' 



Hi = mi(3^^ + p^aW + /3«M(r« - r^°^). 



(25) 



(26) 



Moreover, we have taken in M{x,x') the limit of small Tg and the set of contours in Dk{C) passing from the point 
r^*^ to some (arbitrary) point r^*^). As in [27] one can now integrate over E2,E3 to obtain finally 



G^^'^\E,ri,r'^ J Dk 



(G)(er)(e'r')- 



(27) 



[E — Hi — H2 — H3) 

Prom (27) one obtains equation for the 3q wave function similar to that of qq system, 

{Hi +H2 + H3- E)i){rx, r2, ra) = 0, (28) 
where r^^^ is to be taken at the Torricelli point. In the nonrelativistic approximation mi ^/a one has 



E 



2m, 



+ cr r 



W _r(0)| 



(29) 



IV. PERTURBATIVE CORRECTIONS TO FACTORIZED SOLUTIONS 



The effective Lagrangian (8) and the effective mass operator M{x,y), Eq. (13), do not take into account the 
pcrturbativc interaction between the quarks in the baryon. To this end we separate the gluonic field into a 
background i?^ and perturbative parts, Afj, = Bf^ + and use the 't Hooft identity to integrate in the partition 
function independently over both parts of as was done in [30]. 

We shall use the following representation of gauge transformations 

B^^U+{B^+'-d^)U, a^^U+a^U (30) 

and keep for the background gauge condition [29] , [30] 

D^,{B)a^ = 0, D^{B) =d^- igB^ (31) 

As a result of the perturbative gluon exchange between different quarks in the baryon there will appear an additional 
vertex in the effective Lagrangian [19] 

AL = ^+{x)^f^{x) J' ^+{y)j^My)Mx)a4y))dxdy. (32) 
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In what follows we shall be interested only in the color Coulomb interaction which results from (32) assuming the 
simplest form of gluon propagator and neglecting at first for simplicity the influence of the background field on it, 
namely 

/ {a,ix)aMdix, - ^4) = %^ / ^l:'-''^ = (33) 
J An^ J [x- yy + [Xi - y^Y 47r|x - y| 

Now taking the background into account, one arrives at the picture of the gluon propagating inside the film - the 
world sheet of the string, created by the background between three quark worldlines and the string junction, as is 
shown in Fig. 1. Depending on the choice of r^'^^ we will get in general an effective interaction of a 2-body or 3-body 
nature. Due to the presence of the QCD background the strength of the resulting Coulomb interaction is expected to 
be different from the perturbative OGE contribution and as a result different from the interaction used for example 
in the Breit equation [31]. 

Due to its attractive nature the color Coulomb contribution leads to smaller baryon masses and giving rise to 
composite systems with a smaller radius. As a result the magnetic moments become smaller. In the remaining part 
of the paper we neglect the effect from the Coulomb interaction. To study this a more involved analysis is needed, 
where also the hyperfine interaction has to be included. 




V. BARYON MAGNETIC MOMENTS WITHOUT QUARK CONSTITUENT MASSES 

Since the calculation of magnetic moments as well as baryon masses docs not involve large momentum transfer, one 
can use for that purpose the Hamiltonian equation (28). According to the results of section 4, Hi can be represented 
as 

Hi = m,/3(*) + p(^)a« + (}^^ M^^ {r^'^ - Fq), (34) 

The baryon solution of (34) can be represented as 

*7M = Tf^{hf2h)ea,M].{^^'^ - r(°))V,^|(r(^) - r(o))^/3 (^(3) _ rW), (35) 

where a, h, c and a, /?, 7 refer to color and Lorentz indices respectively and is the flavour index. In what follows we 
shall use only the lowest orbitals (lowest eigenvalues solutions) for quarks and therefore the orbital excitation indices 
are everywhere omitted. The orbital wave function can be decomposed in the standard way 
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^"^^> p \ lF{p)a,VM )"\ lI{p)^hvM ) 



p = r - r(°) (36) 



and the color index is omitted, since the orbital satisfies a "white" (vacuum averaged) equation 

^^iV'4=£^^^V (37) 

Therefore the only remnant of color is the requirement that jm be symmetric in all coordinates besides color. Prom 
Eq. (28) we see that the mass of the baryon, corresponding to Eq. (35), is given by 

3 

To define the magnetic moment one may introduce an external e.m. field A, p(*) p^*^ — eg*^A, A = ^(H x r), and 
calculate perturbatively the energy shift, 

l^E = -/LiH. (39) 

Due to the symmetry of the problem, it is enough to consider only the perturbation of one orbital, say for the first 
quark. 



Hence, denoting ^^^^ = ^ 



H^^Hi+ AHi, AHi = -eWa(i)A. (40) 

(1) 



(A/fr) = -e«(^«*,x(^)*)(^a)Ao 



(1) 



= -e[^){^M*cr^^)Ax^'> + x*^^^ a^^'^ A^^^^). (41) 

Using Eq. (36) and a simple derivation given in Appendix A one obtains for the contribution of the first quark to the 
magnetic moment operator in spin space 



^ 3 



■ j g*{r)fir)rd^rQ*i^a^'^njiM. (42) 

For the lowest orbital J = 5, ^ = 0, M = i, cr ^ a^, one obtains 

/z, = 3/zW = -2eWa(i) J g*{r)f{r)rr^dr, (43) 

where the superscript 1 denotes the contribution of the first quark to the magnetic moment. The normalization 
condition is 



(|5|^ + |/|Vrfr = l. (44) 

Note that everywhere we put r^^) — r(°) = r. In the case of a local linear confining interaction using the Dirac equation 
one can express jj,^^^ through g{r) only (see Appendix A for details) 

^elalf- |,|V(2ar + 3.) 

Constructing the fully symmetrical 3q wave function for the nucleon with total spin up one has for proton 

Kmm = N' {l[u+{l)d.{2) + d.{l)u+{2)]u+{3)- 
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~[d+{l)u-{2) + «_(l)d+(2)]«+(3) - ^[«+(l)«_(2) + «_(l)n+(2)]rf+(3)- 

^[«+(l)d+(2) +rf+(l)«+(2)]u_(3) + |«+(l)u+(2)d_(3)| , (46) 

where = ^ind subscripts (±) refer to the spin projection. In a similar way for the neutron one replaces u d 
and obtains 

= iV' |^[d+(l)«_(2) + ^._(l)d+(2)]d+(3)- 

- i K (2) + rf_ {2)]d+ (3) - 1 (2) + rf_ (2)]u+ (3)- 

i[d+(l)«+(2) + «+(l)d+(2)]rf_(3) + ^rf+(l)d+(2)w_(3)| . (47) 
The matrix elements are computed easily 

o 

where e is the charge of the proton. From Eqs. (48)-(49) one immediately gets the famous relation 
Writing for identical orbitals the magnetic moment as a product 

Mb = 3(*sj/mm|e(^V(l) \^svmm)\ (51) 

where 

X=-^Jg*{r)f{rydr. (52) 

It is clear that inclusion of higher orbitals will change the magnetic moment of proton and neutron, similarly to the 
case of tritium and ^He, where the admixture of the orbital momentum L = 2 changes the magnetic moment by 
7-8%. In our case the orbital momentum is brought by all 3 quarks symmetrically, and these components appear in 
the wave function due to mixing through the tensor and spin-orbit forces between quarks. 

Eqs. (51)-(52) can readily be generalized when the quarks have different orbital wavefunctions. For the single quark 
orbitals we have taken the solution of the Dyson- Schwinger-Dirac equation with nonlocal kernel from Refs. [19,20]. 
Assuming for the field correlator a Gaussian form 

D{u) = D{0)exp{-uy4T^), D{0) = ^ (53) 

with Tg = 0.24 fm the ground state orbital solution is determined. In Table I are shown the calculated ground 
state energy of the orbitals for various flavour states. For the current masses we have used m„ = rud = 5 MeV and 
m, = 200 MeV. 
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Using these orbitals we calculate the nucleon magnetic moment for various values of the string tension a. The 
results are also shown in Table I. From the table we sec that the predictions depends sensitively on the string tension 
<T. Increasing the value of o leads to a larger ground state energy of the orbitals and smaller size of the magnetic 
moment. This in accordance with an analysis, where the small component of the orbital is treated perturbatively. 
Similarly the presence of a Coulomb interaction yields a lower ground state energy of the orbital, resulting in a larger 
value in magnitude of the magnetic moment. Close agreement with the experimental values of the magnetic moment 
is found when a = 0.09 GeV^ . In this case the mass of the nucleon is predicted to be 891 MeV . It is gratifying to 
see, that the magnetic moments are reasonable in the regime where also the predicted mass of the nucleon is close to 
the empirical value. 

TABLE I. Ground state energy eo of the orbitals and the predicted magnetic moments of the nucleons in units of nuclear 
magneton for various values of a. The experimental values are also listed. 



0- (GeV") €o(w,d)(Mey) 


eo(s) (MeV) 


f^proton 


f-^neutron 


0.09 


297 


439 


2.81 


-1.87 


0.12 


342 


482 


2.44 


-1.63 


0.15 


380 


519 


2.20 


-1.46 


experiment 


2.79 


-1.91 


TABLE IL The magnetic moment of the baryons 


in units of nuclear magneton for various values of a. 


Calculations and 


experimental results. 










B 


Mb 




Hb 


Hb 


exp 




a = 0.09 GeV^ 


fT = 0.12 GeV'^ 


fT = 0.15 GeV^ 




P 


2.81 




2.44 


2.20 


2.79 


n 


-1.87 




-1.63 


-1.46 


-1.91 


E" 


-1.03 




-0.89 


-0.79 


-1.16 


E° 


0.85 




0.74 


0.67 




E+ 


2.72 




2.37 


2.14 


2.46 


A 


-0.66 




-0.60 


-0.56 


-0.61 




-0.57 




-0.53 


-0.50 


-0.65 




-1.51 




-1.34 


-1.23 


-1.25 


A++ 


5.62 




4.89 


4.39 


4.52 


A+ 


2.81 




2.44 


2.20 




A° 


0.00 




0.00 


0.00 




A" 


-2.81 




-2.44 


-2.20 




E+* 


3.09 




2.66 


2.37 




E"* 


0.27 




0.21 


0.18 




E"* 


-2.54 




-2.23 


-2.02 






0.55 




0.43 


0.35 






-2.26 




-2.02 


-1.84 






-1.99 




-1.80 


-1.G7 


-2.02 
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The explicit form of "^symm for other baryons arc given in Appendix B. Note, that due to the strange qiiark mass 
their orbitals are different from those of u,d quarks and therefore the decomposition (51) has to be modified. Some 
useful formulas can be found in Appendix B. 

The resulting values for baryon magnetic moments are given in Table II, where they are compared with experimental 
values. Considering the case of cr = 0.12 GeV"^ we see, that there is a rather close agreement with the experimental 
magnetic moments, with the largest deviations found for the nucleon and S~. As discussed for the case of the nucleon 
improvement of the predicted mass of the composite system also leads to magnetic moments closer to the experimental 
values. This applies also for the case of the A-isobar. Hence we may hope that the inclusion of the Coulomb and 
hyperfine splitting interaction will improve the predictions. Moreover, pionic effects are expected to be present. As 
a result significant mesonic current contributions to the magnetic moments may occur. In the next section we study 
the dominant corrections from the pion to the one- and two-body current. 



VI. MESONIC CONTRIBUTIONS 



In this section we carry out in our single orbital model an estimate of the magnitude of the pionic-current corrections 
to the magnetic moment of the nucleon. Due to the quark- coupling to effective mesonic degrees of freedom, one and 
two-body current contributions to the magnetic moments of the baryons exist from the virtual excitations of mesons. 
Assuming as in Ref. [32] that there exists an effective one meson exchange between quarks in the three-quark system 
this leads to meson exchange current contributions to the magnetic moment. The leading correction is due to the 
pion-in-flight and pair term, see Ref. [32]. Effects from the heavier mesons like the p are in general less important. 

Our starting point is the e.m. current matrix element 

M^ = {^\J^{Q)\^), (54) 

where ^ is the 3-quark wavefunction and Q is the photon momentum. 

We first consider the single quark current contribution. For the single quark current operator we use 



^ 3 j7««(i) = 3e(i)^(i) ^(") (55) 

n=2 



and for the wavefunction normalization Eq. (44) for the single particle orbitals is taken. This choice has the nice 
property that the zeroth- component of the current at Q = is found to give the proper charge of the 3-quark system, 
i.e. 

3 

Mo = (vl/|Jo(Q = 0)|vl/) = ^e("). (56) 

n=l 

The result for the magnetic moment, obtained in the previous section can readily be recovered from our single quark 

current matrix element. Following Ref. [33], the magnetic moment can be calculated by taking the curl of the space 
component of the current matrix element in the Breit system. In doing so, the magnetic moment can be deduced 
from the e.m. current as 

= ^^'"'^'^(^ = 0) = - ^ [Vq X M]^ (g = 0), (57) 

where Mp is the proton mass, e the proton charge and Gmag is the Sachs e.m. magnetic form factor. The matrix 
element (57) can easily be evaluated in momentum space. Introducing the Fourier transform of the wavefunction of 
the single quark orbital 

^^(k) = { S'" )=^-l{ ^'^-v'^a^^^o ) P'dp, (58) 
\f{k)Q.jVM J J \t{-ty jv{kp)fip)fljVM J 

with ji the spherical Bessel functions, we may after some algebra reduce Eq. (57) in momentum space to 



H;, = 3/i(^) = 3 (ipsymm e^^V^^) tl^symm) A, (59) 



We thus find, 
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3 

A = 2^ / / d'pd'l n + 1/(^01') 

n=2 

where iV is the normahzation factor 

3 

N= j j d'pd'q n (|5(fc„)P + |/(fc„)|^) . 

n=l 

The momenta are expressed in terms of the Jacobi coordinates as, 



(60) 



Q2=o 



(61) 



k2 



-p+^q+|P 



K = p' + ;^q' + |P', 



(62) 



3' ' 



-P' + 7S^' + 



Imposing the Breit system, P + P' = 0, and momentum conservation gives P' = 
Q. 

Use has been made of the identity 

^jiM{ki) (ki)j(ki)j fljiM{ki)'^ = -5ij. 



Q/2, p' = p and V3(q-q') 



(63) 



with I = and Eqs. (48-49). The magnetic moment expression (51) from the previous section is readily recovered 
when we replace the integration over the Jacobi momenta in Eqs. (60-61) by Y\n=i dkn- 

We now turn to the pionic two-body current contributions, assuming a 75 theory. The resulting pion-in-flight and 
pair current operators, shown in Fig. 2 are given respectively by 



t(23) _ _2,V„2 ^(2) (3) / (2) (3)\ ^ A 

75 (T XT ((A-iQ)Vm=)((A+iQ)Vm= 

((A-iQ)2+A2)((A-HiQ)2+A2) V ^ (A+i(4)2+Aj ^ (A-iQ)2+A2 ; 



(64) 



and 



.(t^-iT 



4m<, 



7 



(2) 



((A-iQ)%™j) ((A-iQ)2+Aj) 

At 



((A+iQ)^+,r.^) ((A+iQ) = +Aj)^ 



(65) 



In Eqs. (64)-(65) Q is the photon momentum, A = p — p'. A monopole form factor with cutoff mass A^r = 675 MeV 
has been used. The last two terms in the last factor in Eq. (64) correspond to contact terms, which are needed to 
satisfy current conservation. The quark-propagator in Eq. (65) has been replaced by its negative energy part. 



P7 — m + ^p2 -|- 77127° 



4m 



(7° - 1) , 



(66) 



as the positive energy part has already been included in the single quark current matrixelement [35]. Moreover, the 
pair contribution (65) consists of 4 terms where the photon can interact with quark 2 and 3 prior and after the 
pion-quark interaction. 

The photo-pion vertex is described by an effective interaction Lagrangian 



(67) 
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k, 



k' 



k'. 




k' 




(a) 



(b) 



FIG. 2. The diagrams corresponding to the pionic contributions to the current: (a) the pion-in-flight diagram, (b) the pair 
term. The bound state of the quarks is represented by the blobs at the beginning and the end of the diagrams. 

Prom the 2-body operators J26, (64-65) we may write down the current matrix element between the 3-quark state 



M26 = 3M^^ 



(1)t(23), 



(68) 



Taking the curl of Eq. (68) the magnetic moment can be determined. The resulting expressions are given in Appendix 

C. As a check using the obtained magnetic moment operators we have determined the exchange magnetic moment 
contribution to the trinucleon system. Our results agree with those obtained by Kloet and Tjon [33]. 



TABLE III. The single quark current contribution /i^^ to the magnetic moment in units of nuclear magneton, together with 



the two-body corrections and the anomalous correction 5/x]v arising from the pion one-loop diagrams. Also are shown the total 
combined prediction of our calculations and the experimental results 





N 




(mr-y) 


(NN-y) 

Mat 




tot 
fJ-N 


exp 








a = 0.09 GeV 








P 


2.81 


0.20 


-0.21 


0.12 


2.92 


2.79 


n 


-1.87 


-0.20 


0.21 


-0.16 


-2.02 


-1.91 








a = 0.12 GeV 








P 


2.44 


0.19 


-0.18 


0.11 


2.56 


2.79 


n 


-1.63 


-0.19 


0.18 


-0.14 


-1.78 


-1.91 








a = 0.15 GeV'-' 








P 


2.20 


0.18 


-0.16 


0.10 


2.32 


2.79 


n 


-1.46 


-0.18 


0.16 


-0.13 


-1.61 


-1.91 
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To get an estimate of the exchange current contributions in the 3-quark case wc have used for the couphngs and 
cut off mass the values from Ref. [32]. They are taken to be g^^^/An = 0.67. The results for the magnetic moments 
are shown in Table III. Our estimates are in strong disagreement with those obtained in Ref. [32]. The pion- in-flight 
contribution is substantially smaller than found in Ref. [32] using the chiral constituent model [34]. This may be 
partially due to the 3-quark wavefunction used, which has a matter radius smaller than in our case. Moreover, it 
contains only nonrclativistic components. The pair contribution is found to be comparable to the pion-in-flight term, 
leading to an almost cancellation of the mesonic current pionic contributions. 





(a) 



(b) 



FIG. 3. The diagrams contributing to the anomalous magnetic moment of the single quark 

The presence of mesonic degrees of freedom will modify the single quark current. The resulting e.m. current 
operator can in general be characterized by a large number of off-shell form factors [36]- [38], which reduces to 2 when 
we assume that the initial and final quark is on-mass shell. Using this approximation we may estimate the resulting 
anomalous magnetic k term due to the the pionic contributions. Near = we have 



Jiqq _ eg 7^ + ^^2M ^^^^^ 



(69) 



where Kg = Ks + k^Tz for the u,d-quark. The k coefficients can be determined in a simple model, assuming that the 
loop corrections are given by only the one-loop pionic contributions to the e.m. vertex. Similarly as in the two-body 
current case we approximate the single quark orbital by a free quark propagation with a constituent mass given by 
the ground state orbital energy. With the above simplifying assumptions the calculation amounts to calculating the 
magnetic moment contributions of the diagrams shown in Fig. 3. Using the same cutoff mass regularization as for the 
two-body currents we find for the anomalous magnetic moment in units of the nuclear magneton, 



. 2 4Mp 



4 {p ■ kY - p'^P 



{2TTr [fc2 - 2pk + ie] [fc2 



A2 



k^-Al 



1 



(70) 



and 



1 - T, 2Mp 
' 2 3m3 



A{p-kY -p'^k'^ 



{2t:Y [k^ - 2pk + ief [k^ - ml + ie] \k^ - Al 



(71) 



where p is the momentum of the quark. For details we refer to Appendix D. Eq. (70) corresponds to the coupling of 
the photon to the pion, Eq. (71) to the coupling of the photon to the quark. 
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TABLE IV. The quark anomalous magnetic moments in units of nucleon magneton in the one-loop approximation for various 
string tension a. The first set is the prediction for only the pion loops, while the second set is with both pion and kaon loops 
included. 



a {GeV^) 




Kd 






pion loops 






0.09 


0.101 


-0.160 


0.0 


0.12 


0.092 


-0.140 


0.0 


0.15 


0.085 


-0.126 


0.0 




pion and kaon loops 






0.09 


0.132 


-0.151 


-0.034 


0.12 


0.121 


-0.133 


-0.032 


0.15 


0.112 


-0.120 


-0.031 



In Table IV are shown the calculated anomalous magnetic moments of the u, d and s quarks for A = 675 MeV for 
various choices of a. Clearly, the results depend on the constituent quark masses. These are given in Table I for the 
considered string tensions. 

Using Eq. (57) the K-term in Eq. (69) yields a nucleon magnetic moment correction 

6lJ.z = = 3 {tpsymm \l^g{'i-)(^z{'i-)\ i'symm) Aq (72) 

with 

Jr^dr{\g\^ + \f\^) ^'''> 

In Table III the predictions for the nucleon are shown including also the one-pion loop contributions (72) and two-body 
currents. Our results obtained for the one-loop corrections are smaller than reported by Glozman and Riska [39]. This 
is due to the inclusion of the lower component in the single quark orbitals. Neglecting these we recover the results 
of Rcf. [39] . From Table III we sec that the proton and neutron magnetic moment is in reasonable agreement with 
experiment for a string tension of cr = 0.1 GeV^. For this value of the string tension the model predicts a nucleon 
mass of 940MeV, remarkably close to the empirical value. The anomalous magnetic moment contributions are found 
to be of the order of 10%. 

Due to the one-loop contributions the magnetic moments of the other baryons are modified. Corrections from kaon 
loops have also been considered. Because of the larger kaon mass the contributions are expected in general to be 
smaller in size than those of the pion loops. In Table IV the calculated anomalous moment of the strange quark due 
to the kaon one-loop corrections are given. In the calculations a cutoff mass of A = 675 MeV has been used. The 
isoscalar and isovector anomalous magnetic moment pieces are also changed by the kaon loop contributions. Prom 
Table IV we see that the kaon loop contributions are indeed smaller in magnitude as compared to the pion loop ones. 
The full results for the magnetic moments of the baryon octet and decuplet, including the pionic exchange currents 
and the pion and kaon one- loop contributions are summarized in Table V. For the value of the string tension cr = 0.1 
the overall agreement with the experimental data is reasonable. From the table we see that the anomalous magnetic 
moment contribution leads to an improvement of the predictions. 
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TABLE V. The magnetic moment ^.b of the baryon octet and decuplet in units of nuclear magneton, including the anomalous 
contribution 6fj,B arising from the pion and kaon one-loop diagrams and the pion exchange corrections for different string tension 
a. Also are shown the experimental results 



B 


a = 0.09 GeV^ 


Ms 


SfXB 

0- = 0.12 GeV"^ 


fJ.B 


SfJ,B 

a = 0.15 GeV'^ 


jJ-B 


exp 


P 


0.15 


2.95 


0.14 


2.59 


0.12 


2.34 


2.79 


n 


-0.16 


-2.02 


-0.14 


-1.78 


-0.13 


-1.61 


-1.91 




-0.13 


-1.16 


-0.11 


-1.00 


-0.10 


-0.89 


-1.16 


S° 


0.00 


0.85 


0.00 


0.74 


0.00 


0.67 




E+ 


0.12 


2.84 


0.11 


2.48 


0.11 


2.25 


2.46 


A 


-0.02 


-0.68 


-0.02 


-0.62 


-0.02 


-0.58 


-0.61 


c— 


0.00 


-0.57 


0.00 


-0.53 


0.00 


-0.50 


-0.65 




-0.06 


-1.57 


-0.05 


-1.39 


-0.05 


-1.28 


-1.25 


A++ 


0.26 


5.88 


0.24 


5.13 


0.22 


4.61 


4.52 


A+ 


0.07 


2.88 


0.07 


2.51 


0.07 


2.27 




A" 


-0.11 


-0.11 


-0.10 


-0.10 


-0.08 


-0.08 




A" 


-0.30 


-3.11 


-0.26 


-2.70 


-0.24 


-2.44 




E+* 


0.15 


3.24 


0.14 


2.80 


0.13 


2.50 




E°* 


-0.04 


0.23 


-0.03 


0.18 


-0.03 


0.15 




E-* 


-0.22 


-2.76 


-0.20 


-2.43 


-0.18 


-2.20 






0.04 


0.59 


0.04 


0.47 


0.03 


0.38 






-0.14 


-2.40 


-0.13 


-2.15 


-0.12 


-1.96 




n- 


-0.07 


-2.06 


-0.06 


-1.86 


-0.06 


-1.73 


-2.02 
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VII. CONCLUSION 

We have written down the general effective quark Lagrangian as obtained from the standard QCD Lagrangian by 
integrating out the ghionic degrees of freedom. Considering the baryon Green's function, neglecting gluon and meson 
exchanges, we find in lowest order of the approximation scheme that it is given by a product of 3 independent single 
quark Green's functions. As a result the Hamiltonian can be written as a sum of three quark terms, where the single 
quark solutions satisfy the Dyson-Schwinger equation with a nonlocal kernel. 

The nonlinear equation for the single quark propagator S (attached to the string in a gauge-invariant way) has been 
solved in the Gaussian correlator approximation. The resulting 3-quark wavcfunction has been used to determine the 
magnetic moments of the baryons. This has been done for both the octet and decuplet of the SU(3) flavour group. 

Comparing the predictions we find that the magnetic moments are mostly in close overall agreement with the 
experiment for a string tension of o" = 0.1 GeV^. We find, that the predicted magnetic moment of the nucleon is 
improved substantially once we choose a string tension to give a reasonable nucleon mass. The same applies for the 
A-isobar. Eff'ects due to the presence of virtual mesons are in general expected to be important. We have estimated 
the pionic one- loop and one pion exchange contributions to the magnetic moment. The single quark corrections from 
pionic loops are found to be of the order of 10%, whereas the total effect of 2-body current contributions are predicted 
to be small, to be contrasted to the results of Ref. [32]. This is due to the cancellation of the pion in flight and 
pair term in the present model. Because of the anomalous magnetic contributions there seems to be somewhat an 
improvement of the predictions. 

Our results for predictions of the magnetic moments of baryons are encouraging, but are in need of including higher 
order corrections. In particular, the mass spectrum obtained from our lowest order approximation does not contain 
the N — A mass splitting. This is due to neglecting contributions like the hyperflne interaction arising from the one 
gluon interaction. It is clearly of interest to investigate how the magnetic moments are changed when effects from 
color Coulomb and hyperflne interaction are accounted for. 
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APPENDIX A: MAGNETIC MOMENT CALCULATION IN COORDINATE SPACE 

The one-quark contribution to the magnetic moment can be written as in (45) 

{AH,) = -e« J (^(1) V(i)Ax(^^ + x^'>a^'^A^^'^)d\ (Al) 

where A = |(H x r) is the vector potential of external constant magnetic field. 

Inserting in (A.l) ip^^^ = g{r)fljiM and x^^^ = ifi'>')^jl'M, and taking into account that flji'M = — (fn)f2j;M, one 
easily obtains 

(AH,) = -ie« I dMg*f + rg)rn*i^{ian){ull) - o-H}0,,m (A2) 

Eq.(A.2) contains the matrix clement J dnQ*ijyjnink^j[M , which simplifies when Z = 0, so that {niUk) = ^Sik- 
In this case one obtains, taking into account relation {AHi) = AE = — /Lt'^^-'H, 

H^'^ = -^eW J {g*{r)f{r) + r{r)g{r))rd'rn*iM<y^'^njiM = 

= _^(i)^eW J Re{g*{r)f{r)ydr (A3) 

In the case of a local scalar potential U{r) one can further express /(r) through g{r) using the Dirac equation for the 
one-quark state 

^^(^^ = . + mU(r) (|:(^^^ + >) 
Introducing (A. 4) into (A. 3) and integrating by parts one obtains 

For U{r) = ar one obtains Eq.(49). 
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APPENDIX B: MAGNETIC MOMENT OF THE MULTIPLET 



In this Appendix the calculation of the nucleon magnetic moment is generalized to the baryon octet and decuplet. 

By analogy with the fully symmetrical 3q wave function for the nucleon, Eqs. (46-47), wave functions for the baryon 
multiplets can be formulated. The flavor octet with total spin 1/2 up becomes, 



XllP 

symm 



symm 



V2 



{2d-U+u+ — U-d+u+ — d+U-U+ + 2u+d-U+ 



—u+U-d+ — U-U+d+ — u+d+U- — d+u+U- + 2u+u+d-} , 

{2?i_d+(i+ — d-U+d+ — u+d-d+ + 2rf+M_rf+ 

6 



E+ 

symm 



d^d^u^ — d+u+d- — u^d^d^ + 2d^d^U-} . 



symm 



^ symm 



— d^d-U^ 

V2 
6 

—U+U-S+ — — U+S+U^ — S^U+U- + 2u+U+S-} , 

— {M+d_s+ + + s+(i_M+ + s+M_(i+ — 2u_|_s_d+ — 2d+S-U^ 
6 

+u_(i+s+ + (i_u+s-|_ — 2s_(i-)_u+ — 2s^u+d+ + u_s+c?+ + c?_s+u_|_ 
—2u+d+s- — 2d+u+s- + s+d+u- + s+u+c?- + u+s+d- + d+s+u_} , 
V2 



symm 



symm 



—— {2s_d+(i+ — (i_s+ci+ — s+(i_(i+ + 2(i+s_ci+ 
6 

— d+d_s+ — d_d+s+ — — s+d+d_ + 2d+d+s_} , 

6 

— s+rf_u+ + s+u_rf+ + s+d+u_ — s+u+rf- — u+s+rf_ + rf+s+u_} , 
\/2 

-— {2m_S+S+ — S_M+S+ — U+S_S+ + 2s_|_'U_S_|_ 

6 

— S+S_U+ — S_S+'U+ — S+M+S- — U+S+S_ + 2s+S+'U_} , 

x/2 



{2(i_s+s^ 



where the subscripts (±) refer to the spin projection. For the flavor decuplet with total spin 3/2 up we have 

T A + + 



A+ 

symm 

A° 

symm 

A- 
symm 



symm 



V3 
1 

1 

1 
1 
1 
1 

7! 

S+S+S+. 



{u-^-u+d-^- + M+rf+u+ + , 



+ + u-^-d+d-^-} , 

{d+d+s+ + d+s+d+ + s+d+d+} , 
{s+s+rf+ + s+d+s+ + d+s+s+} , 



(Bl) 

(B2) 
(B3) 
(B4) 

(B5) 
(B6) 
(B7) 

(B8) 

(B9) 
(BIO) 

(Bll) 

(B12) 
(B13) 

(B14) 

(B15) 

(B16) 

(B17) 
(B18) 
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These fully symmetrical wave functions Eqs. (B1)-(B18) can be written symbolically as, 

^^M,symm = ^fjifMWa'i'i'^^'- (B19) 

As the orbital of the s-quark is heavier than the u- and d-quark orbitals Eq. (51) has to be split up in contributions 
from the u, d-quark and from the s-quark. Using the symmetrical wavefunction Eq. (B19) this is realized by writing, 

M.=3m(i)=3 ^ {rf^'{fif2f3)^l^^^'^^^\e,(l)a,{l)\rf^{^^^^^ (B20) 
/1/2/3 

with, 

H=-\j9Ur)fu{r)r''dr. (B21) 

The flavor index /j can take the values u, d or s. Note that A„ = Ad as the same orbital is taken for the u- and 
rf-quark. Evaluating Eq. (B20) for the difl^erent baryon wave functions Eqs. (B1-B18) results in the expressions in 
Table VI. 



TABLE VI. The matrix elements of the e.m. current for the baryons. 



X 

P 
n 

E+ 

S° 

E" 

A 
=0 



A+ 

A" 
A- 



27/ 



/'-V/'-> 

T\ 

'2 \ 

^ \ \ 
27 -r- 27^^ 

2 X ^ ±\ 

_"4 \ I 1 \ 
27 '^u "T 27 

-±\ - ^\ 

— ■* X -I- 1 X 
27 '^s 27^" 

2a 

q Alt 



- 1 \ 



ix - 

q Au 

l^ 

q 

2X 

~ Q Au 



1 \ 

g As 

q As 
1 \ 

- gAs 

— q As + cjAu 

-lA. - iA„ 
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APPENDIX C: PIONIC TWO-BODY CONTRIBUTION TO THE MAGNETIC MOMENT 



In this Appendix the pion-in-flight and pair contributions to the magnetic moment of the nucleons are given. 
Following Rcf. [33] these contributions arc determined by taking the curl of the pionic two-body currents Eq. (68). 
The 3-quark state ^' is given by the product of three single quark orbitals Eq. (35). Because of symmetry considerations 
it suffices to calculate the magnetic moment contribution of pion exchange between say the second and the third quark 
only and multiply the result by a factor of 3 to include the contribution of the other possible permutations of quark 
pairs. 

Considering the pion-in-flight contribution first (Eq. (64)), taking the curl gives rather long expressions which can 
be divided into two parts, 

The first part gives the larger contribution and can be written as, 

'^^-^^.^ I ^^.W^^^(|^(^OP + l/>OlO X (C2) 
|5(fc^)/>2)5(fc^)/(fc3)^ (pA - p,A.) 

+9{k'2)f{k2)hk'^y9{ki)-^^ ((2p - p' + V3q)A - (2p, -p^ + V3g,)A,) 
+M)9{k2)gik's)f{k3)-^^ ((2p - p' - \/3q)A - {2p, - p'^ - \/3g,)A,) 
+f{k',)g{k2)m)g{ks)-^ (p'A -p;A,) 



3/^2 ^3 



^^^A^mn/ A^ 



A^ + A^y Va^ + a2 

The second part comes from the curl applied to the wavefunctions 

^2 



{^f(''^ys{k'^)f{ks) (k2 X ks)^ - ^f{k2)mrg{ks) (k2 X t;)^ 

^^9{k2)m)hkz) (k; X ka)^ + ^-^~9{k2)hk'Mkz) x kQ^ 



+f{k'2r9{k2)m)Kk3)y (k2 X ka)^ -~f{k'2)9{k2)f{k':,)g{k^)y (k^ x 



X 1 + 2 



2 

A^ + mn/ A2 



A^ + A2y Va^+a^ 

The normalization factor N is the same as used before in the single quark current contribution (Eq. (61)). The 
momenta are expressed in terms of the Jacobi coordinates Eqs. (62) again, but from imposing the Breit system and 
momentum conservation we now get P' = — P = Q/2 and 2v^(q' — q) = Q. In writing down these expressions use 
has been made of the spin-isospin operator sandwiched between the fully symmetric wavefunctions in spin-isospin and 
orbital space of the 3 quarks. 



\Tsymm V ^ J y^symm / \Tsymm V / 2 * J Tsymm/ 



2^ 



(C4) 



It should be noted that the spin-isospin factor (C4) is identical to that found for the tri-nucleon case. For all the other 

baryon wavefunctions given in Appendix B the matrix element of the considered two-body e.m. operators vanish, 
because of the isospin structure of the e.m. operator. Hence the considered two-body currents contribute only to the 
magnetic moment of the proton and neutron. 
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The second part is a relativistic effect which enlarges the values by about 10% and which vanishes in the static 
limit as is shown at the end of this section. 

In the same way the pair term can be analyzed. We find 



with 



1 fl Ak2 - A^(/C2)z V/;/^^^7 \~n'\~fi \ 
' ' 9{k2)f{k2)g{k^)g{k2,) 



k2 \3 



and 



1 f2 Ak^ ^M)A~g^k'2mk2)hk'Mk^)\(-Jh^ (C6) 



(C7) 



k',\3 A'+ml V--— VA^+A^ 



{^^9{k2) {m)f{ks) (k;k3 - {k'2)M.) - f{k',rg{h) (k'^k, - (kUkh)) 

+ ~9{k's)9iks) (^-^^fik2) (k;k2 - {k'2)M.) + ^§^9{k2) {i^2^2 - ikUkh) 

2 



- f{k'2)~9{k2)g{kz)9{k:i)^ (^2^2 - (^2)^(^2).)} 



In the non-relativistic limit the lower component of the wavefunction can be expressed in the upper component as 

/» = -^5(fc), (C8) 

where nriq is the constituent mass of the quark. As a result, the pionic two-body current contributions Eqs. (C1-C3) 
and Eqs. (C5-C7) can be simplified considerably. We obtain for the pion-in-fiight contribution 



^A*^ - eJ^fe^^ /^'g^W + mr9{k2y9{kmkz) 



(AVmj)^ y- + ^ A^; • V A^ 



(C9) 



For the pair term we get 



A"-A= 1 ^ 



(CIO) 



while and vanish. These expressions agree with the results of Refs. [33] and [35]. 
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APPENDIX D: ANOMALOUS MAGNETIC MOMENT CONTRIBUTIONS FROM PION LOOPS 



Our starting point is the e.m. currents, corresponding to the one-loop diagrams shown in Fig. 3, 



(27r) [(p - ky -ml + ie] [fc^ ~ml + ie] {k + Qf -ml + ie 
Al hi / , e~ml , {k + Qf-ml\ 



fc2 - A2 (fc + Q f _ A2 I {k + Q f -kl A;2 - A2 



(Dl) 



and 



j(6) ^ _ • 2 /■ 75 {i>'- /fc + m^) 7^ (^- /; + TO^) 75 / h-l \ , , 

2 y (2;r)"-[(y-A;)2-m2+ie] [(p-fc)2-m2 + ie] [fc2-m2 +ze] U'-A2; ^ ^ 

Since wc have assumed a finite form factor at the irqq vertex, similar as in the two-body current case, the two additicmal 
terms are needed in the last factor of Eq. (Dl) to satisfy current conservation. ^From these currents the anomalous 
magnetic moment has to be extracted. By applying the Gordon decomposition to the current Eq. (69) near Q2 _ q 
it can be seen that the anomalous magnetic moment k is the term proportional to —i^K^ with = + p'^. To 
isolate this term the currents are rewritten by explicit evaluation of the 7-matrix algebra and taking the limit —> 0. 
Using the approximation that the initial and final quark is on-mass shell we obtain 



7(a)-_2io2 e_^i'/_^!^ ^^1^'^" ( Y (1 I 2 



k"^ — ml 



{2wf [fc2 - 2pk + ie] [A;2 - ml + ief \k'^-KJ V fc2 - A2 
-2<<,eT,7^CW, (D3) 



and 



2 



" 2 ^ } (27r)* [fc2 - 2pk + ief [fc2 - m2 + k\ - Al 

1-r 



wlqqe^YCfl (D4) 



As the tensors C'"' depend only on the initial and final momenta they can be written as, 

C« = K^K, + K^Q, + Q^K^ + A« Q^Q^ + Afg^, , (D5) 

where A^*' are Lorentz invariants. It can readily be seen, that only the first term A^*'' contributes to the magnetic 
moment. Substituting Eq. (D5) in Eqs. (D3-D4) and taking the initial and final quark on-mass shell we find for the 
anomalous magnetic moment corrections 

= 8iMpm,52^gT,4"\ (D6) 

n^^) = -AiM,m,glJ-^A^^\ (D7) 
The Lorentz invariant expression can immediately be determined from the tensor C^*2 • We get 

= ^ {^Ki^K^ - K^g''") C« . (D8) 

Inserting Eq. (D8) in (D6-D7) the expressions (70-71) are obtained. 

The kaon one-loop diagrams can be calculated in similar way. The starting point is the expressions Eqs. (Dl- 

D2) again, where the mass of the pion is replaced by the mass of the kaon and the isospin structure is changed to 
{tz + 3F)/2 and — (| + ^Y) respectively in Eqs. (D1-D2) with Y the hypercharge. The expressions for the anomalous 
magnetic moment due to the kaon loop become. 
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-«5Kqql7-2+^rj-— ^ / -— 4— — — o . . ,2 12- \^ 



and 



(27r)'* [A;2 _ 2pA; + m2 - M2 + [fc^ - m^^ + ief \k'^-A\J \ k^ - Aj^ 



(D9) 



« - *5i<rnn n + o ^ ^IZa / TTTT 777 TT ^ 77:^ r7T777 77 712 12" V^-^^J 

















3m^ J 



««V9 3 J J (^2TTf [k^ - 2pk + ml - + le]^ [k^ - m^j. + le] \k^ - Al 

with Mq the mass of the intermediate quark, mg the mass of the initial and final quark. The coupling constant gKqq 
and the cutoff Ak are taken the same as for the pion loop. 
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An effective quark Lagrangian is derived from first principles through bilocal gluon field corre- 
lators. It is used to write down equations for baryons, containing both perturbative and nonper- 
turbative fields. As a result one obtains magnetic moments of octet and decuplet baryons without 
introduction of constituent quark masses and using only string tension as an input. Magnetic mo- 
ments come out on average in reasonable agreement with experiment, except for nucleons and E~. 
The predictions for the proton and neutron are shown to be in close agreement with the empirical 
values once we choose the string tension such to yield the proper nucleon mass. Pionic corrections 
to the nucleon magnetic moments have been estimated. In particular, the total result of the two- 
body current contributions are found to be small. Inclusion of the anomalous magnetic moment 
contributions from pion and kaon loops leads to an improvement of the predictions. 

I. INTRODUCTION 

The QCD dynamics of qq and 3q systems is governed by two basic phenomena: confinement and chiral symmetry 
breaking (CSB), which should be treated in a fully relativistically covariant way. Confinement is usually introduced 
for static quarks via the area law of the Wilson loop [1] or equivalently through the field correlators in the Field 
Correlator Method (FCM) [2,3]. 

For spinless quarks, or neglecting spin-dependent mass corrections, one can envisage a self-consistent method which 
treats confinement as the area law also for light quarks in a relativistically covariant way. Such method was introduced 
originally in [4] for mesons, in [5] for baryons, and in [6] for heavy-light mesons, and later on in [7] the method was 
generalized taking into account the dynamical degrees of freedom of the QCD string, which naturally appears due to 
the area law. 

As a result Regge trajectories have been found in [7] with the correct string slope {2TTa)~^. It was realized later on, 
that the method used in [4]- [7] can be more generally developed in the framework of the so-called einbein formalism, 
sec [8]- [10]. Spin corrections have been considered in [11] for heavy mesons and in [6] for heavy-light ones. In the 
general case of light quarks spin-dependent correlations have been introduced in [12], and for gluons in [13]. For a 
general review with explicit formulae see [14]. Baryon Regge trajectories have been found in [5]. In all cases the basic 
formalism is the FCM and the Feynman-Schwinger (or world-line) path integral representation [3,15,16] which is well 
suited for relativistic quarks when spin is considered as a perturbation. 

The main difficulty which was always present in the method, was the perturbative treatment of spin degrees of 
freedom (which is incorrect, e.g., for the pion) and absence of spontaneous CSB effects in general [17]. Recently a new 
type of formalism was suggested to treat simultaneously confinement and CSB and a nonlinear equation was derived 
for a light quark in the field of heavy antiquark [18]. This equation derived directly from QCD Lagrangian was found 
to produce linear confinement and CSB for the light quark and the explicit form of the effective quark mass operator 
M(x,y) was defined obeying both these properties. 

The eigenvalues and eigenf unctions of the nonlocal and nonlinear equations have been determined and a nonzero 
condensate was computed in [19], confirming that CSB is really present in the equations. In an additional study [20] it 
was demonstrated that magnetic field correlators do not contribute to the large distance confinement, however strongly 
modify the confinement for lowest levels and heavy-light masses corrected in this way are favourably compared in [20] 
to the experiment and results of other calculations. 

Moreover, it was shown in [21] that lattice data strongly support the dominance of the Gaussian (bilocal) correlator, 
estimating the correction due to higher correlators to 1-2%. Since the method of [18] is quite general and allows to 
treat also multiquark systems, it can be applied to the qq and 3q systems, to find dynamical equations for them, 
which contain confinement and CSB [22]. To make these equations tractable, one systematically exploits the large Nc 
limit, and mostly confine ourselves to the simplest field correlators - the so-called Gaussian approximation; it was in 
particular shown in [18] that the sum over all correlators does not change the qualitative results. However, the kernel 
of equations becomes much more complicated. 

In the present paper we study the baryon magnetic moments based on the derived effective Lagrangian without 
constituent quark masses. The paper is organized as follows. In Section 2 the general effective quark Lagrangian 
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from the standard QCD Lagrangian is obtained by integrating out gluonic degrees of freedom, and the nonhnear 
equation for the single quark propagator S (attached to the string in a gauge-invariant way) is derived, following the 
procedure in [22]. Section 3 is devoted to the baryon Green's function, which can be expressed in the lowest order of 
our approximation scheme (neglecting ghion and pion exchanges) in terms of 3 independent quark Grecm's fmiction, 
resulting in a Hamiltonian as a sum of three quark terms. In section 4 the next order approximation is written down 
when perturbative gluon exchanges are taken into account, including the nonperturbative interaction between quarks 
violating the factorized form of the zcroth order approximation. The next section is devoted to the calculation of 
magnetic moments of baryons both in octet and decuplet representations of SU(3) flavour group. In section 6 we 
discuss the corrections to magnetic moments due to pion exchange contributions. 

II. EFFECTIVE QUARK LAGRANGIAN 

As was discussed in the previous section, one can obtain an effective quark Lagrangian by averaging over background 
gluonic fields. We shall repeat this procedure following [18] now paying special attention to the dependence on the 
contour in the definition of contour gauge, and introducing the operation of averaging over contour manifold. The 
QCD partition function for quarks and gluons can be written as 

Z = J DADi,D7p+exp[Lo + Li+Lint], (1) 

where we are using Euclidean metric and define 

Lo = -\j d'x{F;^)\ (2) 

Lx = -ij ftp+{x){d + mf) ^^{x)d'^x, (3) 

^■ip+{x)gA{x) ^'■ip{x)d^x. (4) 

Here and in what follows ^ipaa denotes quark operator with flavour /, color a and bispinor index a. 

To express An{x) through F/^v one can use the generalized Fock-Schwinger gauge [23] with the contour C{x) from 
the point a; to a;o, which can also lie at infinity, 

M^) = J/^0i^)^^^ds. (5) 

Now one can integrate out the gluonic field A^(a;), and introduce an arbitrary integration over the set of contours 
C{x) with the weight Dk,{C), since Z is gauge invariant it does not depend on the contour C{x). One obtains 

Z = j DK{C)Di;DtP+exp{Li + L^s}, (6) 
where the effective quark Lagrangian Left is defined as 

expLe. = (exp/ Vi^M.. (7) 



Using the cluster expansion, Leff can be written as an infinite sum containing averages {{A)'') a- At this point one 
can exploit the Gaussian approximation, neglecting all correlators {{A)'') of degree higher than k = 2. Numerical 
accuracy of this approximation was discussed and tested in [21]. One expects that for static quarks corrections to 
Gaussian approximation amount to less than 2-3%. 
The resulting effective Lagrangian is quartic in tp, 

= ^ / d^xdVi'Lix) ^Mx) 'i>ty{y) 'My)Ja/^;j5ix, y) + 0{i^% (8) 
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and J^i, is expressed as 



J^lly{x,y) = g"^ [ ^^du^ [ ^^dUe'-^{Fsu,{u)Fs' 
Jc dxij, Jc dyv Nc' 



'{v)). (10) 



Les, Eq. (8), is written in the contour gauge [23]. It can be identically rewritten in the gauge-invariant form if one 

substitutes parallel transporters ^{x,xo),^{y,xo) (identically equal to unity in this gauge) into Eqs. (8) and (10), 
multiplying each ip{x) and "0(2/) respectively and replacing F{u) in Eq. (10) by <J>(a;, m)F(u)$(u, ccq) and similarly for 
F{v). 

After that Loff becomes gauge-invariant, but in general contour -dependent, if one keeps only the quartic term (8), 
and neglects all higher terms. A similar problem occurs in the cluster expansion of Wilson loop, when one keeps only 
lowest correlators, leading to the (erroneous) surface dependence of the result. The situation here is the same as with 
a sum of QCD perturbation series, which depends on the normalization mass /i for any finite number of terms in the 
series. This unphysical dependence is usually treated by fixing ^ at some physically reasonable value fio (of the order 
of the inverse size of the system) . 

The integration over contours Dk{C) in (6) resolves this difficulty in a similar way. Namely, the partition function Z 
formally does not depend on contours (since it is integrated over a set of contours) but depends on the weight Dk{C). 
We choose this weight in such a way, that the contours would generate the string of minimal length between q and q. 
Thus the physical choice of the contour corresponds to the minimization of the meson (baryon) mass over the class of 
strings, in the same way as the choice of /z = /xq corresponds to the minimization of the dropped higher perturbative 
terms. As a practical outcome, we shall keep the integral Dk{C) till the end and finally use it to minimize the string 
between the quarks. 

Till this point we have made only one approximation —neglected all field correlators except the Gaussian one. Recent 
lattice calculations (see Refs. [24,25]) estimate the accuracy of this approximation at the level of few percents. Now 
one must use another approximation, i.e. assume a large Nc expansion and keep the lowest term. As was shown in 
[18] this enables one to replace in (8) the colorless product -^t/Jbix) ^i^^{y) = tr{ ^ip{x)^{x,Xo)^{xo,y) ^'4'~^{y)) by 
the quark Green's function 



V6/3(x) ^6+ (y) ^ SfgN.S/s^ix, y). (11) 



assumes the form 



^eff =-ij d'^d% f^Pt^ix) fMas{x,y) ^My), (12) 
where the quark mass operator is 

^Mas{x,y) = -J^L/(a;,y)(7^ ^ S{x,y)'y^)aS- (13) 
Prom (12) it is evident that ■''5 satisfies 

{-id^-imf) ^S{x,y)-i j f M{x,z)d'^z ^ S{z,y) = 5^'^\x - y). (14) 

Eqs. (13)-(14) have been first derived in [18]. Prom (6) and (12) one should expect that at large Nc the qq and 

Zq dynamics is expressed through the quark mass operator (13), which should contain both confinement and GSB. 
Indeed, the analysis performed in Refs. [18]- [20] reveals that confinement is present in the long distance form of 
M{x,y), when both distances |x|, |y| of light quark from heavy antiquark (placed at x = 0) are large. 

We shall now make several simplifying assumptions, to clarify the structure of M{x,y). First of all we take the 
class of contours C going from any point x = (.T4.X) to the point (x4,0) and then to (—00) along the X4 axis. For 
this class the corresponding gauge was studied in [26]. Secondly, we take the dominant part of J^i, in (13), namely 
J44, which is proportional to the correlator of color-electric fields. This yields a linear confining interaction, while 
the other components Jik, Ji4,J4i,i = 1,2,3 have been neglected, containing magnetic fields and yielding momentum 
dependent corrections. (It is easy to take into account these contributions in a more detailed analysis). 

The correlator {FF) in (10) can be expressed through the scalar correlator D{x), defined as [2], 



trg^ 

Nr. 



{Fc0{u)^u,v)F^s{v)^v,u)) = D{u - v){Sa',60s - SasS/sj) + 0{Di), (15) 



3 



where the correlator Di, not contributing to confinement, is neglected. As a result one has for M [19,20] 

^Mc^^{x,y)= ■fMW/+ ■fM«ai+ /m(4)74 + ^M«7,. (16) 

Here we have defined 

- = (o-°,)- <>^) 

The dominant part of M, -^M^"^ is linearly growing at large |x|, |y| and in the most simple case of Gaussian form of 
D{x), can be written as 

/MW(a=,t/) = 1 e-'^a|^|5-(3)(x-y) (18) 
where Tg is the gluon correlation length, and ^ is a smeared (5-function, which can be represented as [19,20] 

5(3)(x-y)«earp(-^i^)(^)^ 6 ~ 2T,. (19) 
Here again T„ is the gluon correlation length, which enters -D(w) as -D(m) = -D(0)exp(— ^^). We are now in the 
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position to derive the gg, 2>q Green's function, which will be done in the next section. 



III. EQUATIONS FOR THE BARYON GREEN'S FUNCTION 

Equations for the system can be written in the same way as for the qq system. We again shall assume the large 
TVc limit in the sense, that 1/A^c corrections from qq pairs to the quark Green's function and the effective mass can 
be neglected. We now write down the explicit expressions for N^. = 3. 
The initial and final field operators are 

'^inix,y,z) = eabcr"'^^^aa{x,C{x))i;b0{y,C{y))^c'y{z,C{z)) (20) 

with the notations: a, b, c, are color indices, a, /3, 7 are Lorentz bispinor indices and transported quark operators are 

iJaa{x,C{x)) = {^c{x,x)'llJa{x))a- (21) 

The contour C{x) in can be arbitrary , but it is convenient to choose it in the same class of contours that is used 

in Dk{C) and in the generalized Fock Schwingcr gauge [23]. T°'P'^ is the Lorentz spinor tensor securing proper baryon 
quantum numbers. One can also choose other operators, but it does not influence the resulting equations. In Eq. (20) 
we have omitted flavour indices in F and tp{x,C), to be easily restored in final expressions. 

Using now the effective Lagrangian (12) valid at large Nc, we obtain for the 3q Green's function. 

G^^''\x,y,z\x',y',z') = 
1 J Dn{C)D^DxP+^f,^{x', y', 2')*i(x, y, z)eMLi + L,s). (22) 
Integrating out quark degrees of freedom and neglecting the determinant at large Nc, one has 

^(39) = J DK{C){er){e'T'){S{x, x')S{y, y')S{z, z') + perm}, (23) 

where for simplicity color and bispinor indices are suppressed together with parallel transporters in initial and final 
states. 

One can also define unprojected (without F, F') 3q Green's function Gun^ with 3 initial and 3 final bispinor indices 
instead of projected by F, F' quantum numbers of baryon. Assuming that minimization over contours Dk{C) reduces 
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to the single choice of the contours (the single string junction trajectory minimizing the mass of baryon), one can 
write equation for Gur?^ : 

{-idx - imi - iMi){-idy - im2 - iM2){-idz - ims - iM3)G'^^^ = 

X - x')S^^^ {y - y')5^'^^ {z - z') (24) 

with e.g. MiG = / M{x,u)G{u,x')d'^u. One can simphfy the form (23) for G'^"^^ taking into account that M{x,x') 
actually does not depend on Hence the interaction kernel of G^^'^ does not depend on relative energies, as 

in [27] . Similarly to [27,28] one can introduce Fourier transform of G^^'^ in time components and take into account 
energy conservation E = Ei + E2 + E3. One obtains 

G^^''\E,E2,E3)^ j DK{C){eT){e'T') 

1 



(25) 



[E — E2 — E3 — Hi)[E2 — i?2)(-E'3 — H3) ' 
where we have used the notation 

Hi = m,f3^'^ + p(*)a(*) + (3'^''> M {r^''> - r(°)). (26) 

Moreover, we have taken in M{x,x') the limit of small Tg and the set of contours in Dk{C) passing from the point 
r^*^ to some (arbitrary) point r^^'. As in [27] one can now integrate over E2,E3 to obtain finally 

G(3«)(E,r„r;) ^ J Z)«(G)(er)(eTO^^--^-^^— ^. (27) 

Prom (27) one obtains equation for the 3q wave function similar to that of qq system, 

iHi+H2 + H3-E)^iri,r2,r3)=0, (28) 
where r'"^ is to be taken at the TorriceUi point. In the nonrelativistic approximation nii 3> one has 

y [^^^ + (7]rW-rW]l * = s = E-ymi (29) 
r— f 2mi ^-^ 



IV. PERTURBATIVE CORRECTIONS TO FACTORIZED SOLUTIONS 



The effective Lagrangian (8) and the effective mass operator M{x,y), Eq. (13), do not take into account the 
pcrturbativc interaction between the quarks in the baryon. To this end we separate the ghionic field into a 
background and perturbative parts, = + and use the 't Hooft identity to integrate in the partition 
function independently over both parts of ^4^ as was done in [30]. 

We shall use the following representation of gauge transformations 

B^ ^ U+{B^ + '-d^)U, ^ U+a^JJ (30) 

and keep for the background gauge condition [29] , [30] 

D^{B)a^ = 0, D^{B) =d^- igB^ (31) 

As a result of the pcrturbativc gluon exchange between different quarks in the baryon there will appear an additional 
vertex in the effective Lagrangian [19] 

AL = ^+{x)^l^{x) j' ij+{y)^ti,{y){a^{x)a,{y))dxdy. (32) 
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In what follows we shall be interested only in the color Coulomb interaction which results from (32) assuming the 
simplest form of gluon propagator and neglecting at first for simplicity the influence of the background field on it, 
namely 

/ {a,{x)a4y))d{x, - = / , '1:1- = -4^. (33) 

J Att^ J [x- yy + [Xi - ViY 47r|x - y| 

Now taking the background into account, one arrives at the picture of the gluon propagating inside the film - the 
world sheet of the string, created by the background between three quark worldlincs and the string junction, as is 
shown in Fig. 1. Depending on the choice of r(°) we will get in general an effective interaction of a 2-body or 3-body 
nature. Due to the presence of the QCD background the strength of the resulting Coulomb interaction is expected to 
be different from the perturbative OGE contribution and as a result different from the interaction used for example 
in the Breit equation [31]. 

Due to its attractive nature the color Coulomb contribution leads to smaller baryon masses and giving rise to 
composite systems with a smaller radius. As a result the magnetic moments become smaller. In the remaining part 
of the paper we neglect the effect from the Coulomb interaction. To study this a more involved analysis is needed, 
where also the hyperfine interaction has to be included. 




V. BARYON MAGNETIC MOMENTS WITHOUT QUARK CONSTITUENT MASSES 

Since the calculation of magnetic moments as well as baryon masses does not involve large momentum transfer, one 
can use for that purpose the Hamiltonian equation (28). According to the results of section 4, Hi can be represented 

as 

Hi = m^/3(') + p«a« + /?«MW(r(^) - ro), (34) 

The baryon solution of (34) can be represented as 

^JM = Tf^{hf2h)ea,ci^t{v('^ - rW)^,^|(r(2) - r(o))V,/3(r(3) _ ^W), (35) 

where a, b, c and a, /?, 7 refer to color and Lorcntz indices respectively and fi is the flavour index. In what follows we 
shall use only the lowest orbitals (lowest eigenvalues solutions) for quarks and therefore the orbital excitation indices 
are everywhere omitted. The orbital wave function can be decomposed in the standard way 
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and the color index is omitted, since the orbital satisfies a "white" (vacuum averaged) equation 

^^iV'^=£^^^V (37) 

Therefore the only remnant of color is the requirement that ^ jm be symmetric in all coordinates besides color. From 
Eq. (28) we see that the mass of the baryon, corresponding to Eq. (35), is given by 

3 

Mb = Y.^^I (38) 

To define the magnetic moment one may introduce an external e.m. field A, p'^*^ p^*^ — eg ^A, A = ^(H x r), and 
calculate perturbatively the energy shift, 

/^E = -nH. (39) 

Due to the symmetry of the problem, it is enough to consider only the perturbation of one orbital, say for the first 
quark. 



Hi^Hi + AHi, AHi = -e^^^a^^U. (40) 

/ Ji) \ 

Hence, denoting ^^^^ = 



g 

(1) 



(A^.) = -ea)(,a)*,,a>)(X 

= -eW(^(i)v(i)Ax(i)* + x*^'V(i)A^(i)). (41) 

Using Eq. (36) and a simple derivation given in Appendix A one obtains for the contribution of the first quark to the 
magnetic moment operator in spin space 



^ ~ 3 



■ J g*{r)fir)rd^rQ*ij^am^iM. (42) 
For the lowest orbital j = \, Z = 0, M = 5, cr — > cr^, one obtains 

/z, ^ 3m(^) = -2eWa(i) / g*{T)f{r)rr^dr, (43) 



where the superscript 1 denotes the contribution of the first quark to the magnetic moment. The normalization 
condition is 



/(Iffl^ + l/lV 



)r^dr = 1. (44) 

Note that everywhere we put r^^^ — r^*^^ = r. In the case of a local linear confining interaction using the Dirac equation 
one can express /x^'^ through g{r) only (see Appendix A for details) 

^ 4Mi r |g|V(2ar + 3.) 
Constructing the fully symmetrical 3g wave function for the nucleon with total spin up one has for proton 

Kmm = N'[\[u+{l)d.{2) + d.{l)u+{2)]u+{^)- 
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-^[d+(l)«_(2) + «_(l)d+(2)]«+(3) - ^[«+(l)«_(2) + «_(l)«+(2)]d+(3)- 

i (2) + d+{l)u+{2)]u. (3) + (l)«+(2)rf_ (3) I , (46) 

where = and subscripts (±) refer to the spin projection. In a similar way for the neutron one replaces m d 
and obtains 

= Af' |^[rf+(l)«-(2) + «_(l)d+(2)]d+(3)- 

- i [«+ (2) + rf_ (2)]d+ (3) - i (2) + d_ (2)]u+ (3)- 

^[rf+(l)u+(2) +«+(l)d+(2)]rf_(3) + ^d+(l)d+(2)u_(3)| . (47) 
The matrix elements are computed easily 

le^^V^^l*'' ^ = -e f48) 
2 

\ symml^q \^ symml g^' V ^/ 

where e is the charge of the proton. From Eqs. (48)-(49) one immediately gets the famous relation 

^ = (50) 
3 ^ ' 

Writing for identical orbitals the magnetic moment as a product 

HB = H'isymm\e^q^Cr^P\^symm)\ (51) 

where 

X=-\jg*{r)f{r)r^dr. (52) 

It is clear that inclusion of higher orbitals will change the magnetic moment of proton and neutron, similarly to the 
case of tritium and ^He, where the admixture of the orbital momentum L = 2 changes the magnetic moment by 
7-8%. In our case the orbital momentum is brought by all 3 quarks symmetrically, and these components appear in 
the wave function due to mixing through the tensor and spin-orbit forces between quarks. 

Eqs. (51)-(52) can readily be generalized when the quarks have different orbital wavefunctions. For the single quark 
orbitals we have taken the solution of the Dyson-Schwinger-Dirac equation with nonlocal kernel from Refs. [19,20]. 
Assuming for the field correlator a Gaussian form 

D{u) = D{Q)exp{-u-'/AT^^), D{0) = ^ (53) 

with Tg = 0.24 fm the ground state orbital solution is determined. In Table I are shown the calculated ground 
state energy of the orbitals for various flavour states. For the current masses we have used m„ = nid = 5 MeV and 
THs = 200 MeV. 
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Using these orbitals we calculate the micleon magnetic moment for various values of the string tension a. The 
results are also shown in Table I. From the table we see that the predictions depends sensitively on the string tension 
a. Increasing the value of a leads to a larger ground state energy of the orbitals and smaller size of the magnetic 
moment. This in accordance with an analysis, where the small component of the orbital is treated perturbativcly. 
Similarly the presence of a Coulomb interaction yields a lower ground state energy of the orbital, resulting in a larger 
value in magnitude of the magnetic moment. Close agreement with the experimental values of the magnetic moment 
is foimd when a = 0.09 GeV"^. In this case the mass of the nuclcon is predicted to be 891 MeV. It is gratifying to 
see, that the magnetic moments are reasonable in the regime where also the predicted mass of the nucleon is close to 
the empirical value. 

TABLE I. Ground state energy eo of the orbitals and the predicted magnetic moments of the nucleons in units of nuclear 
magneton for various values of a. The experimental values are also listed. 



a (G 


'\"-) cu{u, 


d) (Mc\') 






fiprolu,, 


I'm. ulron 


0.09 


297 


439 




2.81 


-1.87 


0.12 


342 


482 




2.44 


-1.63 


0.15 


380 


519 




2.20 


-1.46 


experiment 




2.79 


-1.91 


TABLE II. The magnetic moment of the baryons in units of nuclear magneton for 


various values of a. 


Calculations and 


experimental results. 












B 


Hb 








Hb 


exp 




a = 0.09 GeV^ 


0- = 0.12 GeV"^ 




0- = 0.15 GeV'^ 




P 


2.81 




2.44 




2.20 


2.79 


n 


-1.87 




-1.63 




-1.46 


-1.91 


S" 


-1.03 




-0.89 




-0.79 


-1.16 


s° 


0.85 




0.74 




U.D( 




E+ 


2.72 




2.37 




2.14 


2.46 


A 


-0.66 




-0.60 




-0.56 


-0.61 




-0.57 




-0.53 




-0.50 


-0.65 




-1.51 




-1.34 




-1.23 


-1.25 


A++ 


5.62 




4.89 




4.39 


4.52 


A+ 


2.81 




2.44 




2.20 




A" 


0.00 




0.00 




0.00 




A- 


-2.81 




-2.44 




-2.20 




S+* 


3.09 




2.66 




2.37 




S°* 


0.27 




0.21 




0.18 






-2.54 




-2.23 




-2.02 






0.55 




0.43 




0.35 






-2.26 




-2.02 




-1.84 




n- 


-1.99 




-1.80 




-1.67 


-2.02 



9 



The explicit form of "^symm for other baryons are given in Appendix B. Note, that due to the strange quark mass 
their orbitals are different from those of u,d quarks and therefore the decomposition (51) has to be modified. Some 
useful formulas can be found in Appendix B. 

The resulting values for baryon magnetic moments are given in Table II, where they arc compared with experimental 
values. Considering the case of cr = 0.12 GeV^ we see, that there is a rather close agreement with the experimental 
magnetic moments, with the largest deviations found for the nucleon and E~. As discussed for the case of the nucleon 
improvement of the predicted mass of the composite system also leads to magnetic moments closer to the experimental 
values. This applies also for the case of the A-isobar. Hence we may hope that the inclusion of the Coulomb and 
hyperfine splitting interaction will improve the predictions. Moreover, pionic effects are expected to be present. As 
a result significant mesonic current contributions to the magnetic moments may occur. In the next section we study 
the dominant corrections from the pion to the one- and two-body current. 



VI. MESONIC CONTRIBUTIONS 



In this section we carry out in our single orbital model an estimate of the magnitude of the pionic-current corrections 
to the magnetic moment of the nucleon. Due to the quark-coupling to effective mesonic degrees of freedom, one and 
two-body current contributions to the magnetic moments of the baryons exist from the virtual excitations of mesons. 
Assuming as in Ref. [32] that there exists an effective one meson exchange between quarks in the three-quark system 
this leads to meson exchange current contributions to the magnetic moment. The leading correction is due to the 
pion-in-fiight and pair term, see Ref. [32]. Effects from the heavier mesons like the p are in general less important. 

Our starting point is the e.m. current matrix element 

M^ = (*|J^(g)|*), (54) 

where ^ is the 3-quark wavefunction and Q is the photon momentum. 

We first consider the single quark current contribution. For the single quark current operator we use 



j;"" ^ 3J;'"'(1) = 36^7^'^ n ^o"^ (55) 

n=2 



and for the wavefunction normalization Eq. (44) for the single particle orbitals is taken. This choice has the nice 
property that the zeroth-component of the current at Q = is found to give the proper charge of the 3-quark system, 
i.e. 



Mo = (*|Jo(Q = 0)|*) = ^4"^ (56) 

n=l 

The result for the magnetic moment, obtained in the previous section can readily be recovered from our single quark 
current matrix element. Following Ref. [33], the magnetic moment can be calculated by taking the curl of the space 
component of the CTirrent matrix element in the Breit system. In doing so, the magnetic moment can be deduced 
from the e.m. current as 

= ^^GmagiQ = 0) = - ^ [Vq X M]^ (Q = 0), (57) 

where Mp is the proton mass, e the proton charge and Gmag is the Sachs e.m. magnetic form factor. The matrix 
element (57) can easily be evaluated in momentum space. Introducing the Fourier transform of the wavefunction of 
the single quark orbital 

\f{k)njVM J J \i{-iy 3v{kp)f{p)VLjVM J 

with ji the spherical Bessel functions, we may after some algebra reduce Eq. (57) in momentum space to 



Hz — ^ — 3 (^tpsymm 

We thus find, 



i>symm)\ (59) 
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^ f /rfv«n(i5(Mp+i/(M 



where iV is the normahzation factor 

N= [ [ SpdPq n {\~9{kn)? + |/(fc„)r) 
The momenta are expressed in terms of the Jacobi coordinates as, 



(60) 



Q2=o 



(61) 



k3 = -p + ^q+iP, 



K = p 
k' =- 



' + 73^' + !?'' 



(62) 



Imposing the Breit system, P + P' = 0, and momentum conservation gives P' = — P = Q/2, p' = p and •\/3(q — q') 

Q 

Use has been made of the identity 



^jimiki) (ki)j(ki)j- VljiMih) 



(63) 



with I = and Eqs. (48-49). The magnetic moment expression (51) from the previous section is readily recovered 
when we replace the integration over the Jacobi momenta in Eqs. (60-61) by rin=i '^^n- 

We now turn to the pionic two-body current contributions, assuming a 75 theory. The resulting pion-in-flight and 
pair current operators, shown in Fig. 2 are given respectively by 



j(23) _ _2,Vo2 ^(2) (3) / (2) (3)\ A 

Jj^n- ^«e5,gq75 75 XT j. ((A-iQ)%™^,) ((A+iQ)^- 



((A-ig)2+A2)((A+ig)2+A2) \} ^ (A+iQ)2+A= ^ (A-ig)2+A2 



(A-iQ)^+. 



(A+iQ)^+. 



(64) 



and 



.(2) 



f^O ,\(3) 



4m„ 



-T 



(7"-l)'^' 
4ni„ 



{(A-^QY+ml) ((A-iQ)2+A|)^ 



7^ 



Ai 



((A+iQ)Vm^) ((A+1Q)^+A2)^ 



(65) 



In Eqs. (64)-(65) Q is the photon momentum, A = p — p'. A monopole form factor with cutoff mass = 675 MeV 
has been used. The last two terms in the last factor in Eq. (64) correspond to contact terms, which are needed to 
satisfy current conservation. The quark-propagator in Eq. (65) has been replaced by its negative energy part. 



2^p2 + m2 pO + ^p2 + ; 



4m 



(7° - 1) 



(66) 



as the positive energy part has already been included in the single quark current matrixelement [35]. Moreover, the 
pair contribution (65) consists of 4 terms where the photon can interact with quark 2 and 3 prior and after the 

pion-quark interaction. 

The photo-pion vertex is described by an effective interaction Lagrangian 



(67) 
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(a) fb) 



FIG. 2. The diagrams corresponding to the pionic contributions to the current: (a) the pion-in-flight diagram, (b) the pair 
term. The bound state of the quarks is represented by the blobs at the beginning and the end of the diagrams. 

From the 2-body operators J 2b, (64-65) we may write down the current matrix element between the 3-quark state 

M2b = 3MW =^jjj J d'pd^qH'^j'i,^^^. (68) 

Talcing the curl of Eq. (68) the magnetic moment can be determined. The resulting expressions are given in Appendix 
C. As a check using the obtained magnetic moment operators we have determined the exchange magnetic moment 
contribution to the trinucleon system. Our results agree with those obtained by Kloet and Tjon [33]. 

TABLE III. The single quark current contribution /it^^ to the magnetic moment in units of nuclear magneton, together with 
the two-body corrections and the anomalous correction arising from the pion one-loop diagrams. Also are shown the total 
combined prediction of our calculations and the experimental results 



N 






(NN-y) 




tot 
f^N 


exp 








a = 0.09 GeV'^ 








P 


2.81 


0.20 


-0.21 


0.12 


2.92 


2.79 


11 


-i.87 


-0.20 


0.21 


-O.K) 


-2.02 


-1.91 








a = 0.12 GeV 








P 


2.44 


0.19 


-0.18 


0.11 


2.56 


2.79 


n 


-1.63 


-0.19 


0.18 


-0.14 


-1.78 


-1.91 








a = 0.15 GeV 








P 


2.20 


0.18 


-0.16 


0.10 


2.32 


2.79 


n 


-1.46 


-0.18 


0.16 


-0.13 


-1.61 


-1.91 
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To get an estimate of the exchange current contributions in the 3-quark case we have used for the couphngs and 
cut off mass the values from Ref. [32]. They are taken to be ggg^/An = 0.67. The results for the magnetic moments 
are shown in Table III. Our estimates are in strong disagreement with those obtained in Ref. [32] . The pion- in-flight 
contribution is substantially smaller than foiuid in Ref. [32] using the chiral constituent model [34]. This may be 
partially due to the 3-quark wavefunction used, which has a matter radius smaller than in our case. Moreover, it 
contains only nonrelativistic components. The pair contribution is found to be comparable to the pion-in-flight term, 
leading to an almost cancellation of the mesonic current pionic contributions. 





(a) 



(b) 



FIG. 3. The diagrams contributing to the anomalous magnetic moment of the single quark 

The presence of mesonic degrees of freedom will modify the single quark current. The resulting e.m. current 

operator can in general be characterized by a large number of off-shell form factors [36]- [38], which reduces to 2 when 
we assume that the initial and final quark is on-mass shell. Using this approximation we may estimate the resulting 
anomalous magnetic k term due to the the pionic contributions. Near = we have 



Jiqq _ eq7p -|- '^9 2M ^^^^^^ 



(69) 



where Kq = Kg + HvTz for the u,d-quark. The k coefficients can be determined in a simple model, assuming that the 
loop corrections are given by only the one-loop pionic contributions to the e.m. vertex. Similarly as in the two-body 
current case we approximate the single quark orbital by a free quark propagation with a constituent mass given by 
the ground state orbital energy. With the above simplifying assumptions the calculation amounts to calculating the 
magnetic moment contributions of the diagrams shown in Fig. 3. Using the same cutoff mass regularization as for the 
two-body currents we find for the anomalous magnetic moment in units of the nuclear magneton. 



4(p- fc) -p2fc2 



3m3 J {2TTf [k^ - 2pk + le] [P - ml + %tf - 



1 + 2 



(70) 



and 



2 1 - 2Mp 



d'^k 



4 {p ■ ky - p^k"^ 



2 3m3 J (27r)'' [k^ - 2pk + ief [k^ - ml + le] \ k 



A2 



(71) 



where p is the momentum of the quark. For details we refer to Appendix D. Eq. (70) corresponds to the coupling of 
the photon to the pion, Eq. (71) to the coupling of the photon to the quark. 
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TABLE IV. The quark anomalous magnetic moments in units of nucleon magneton in the one-loop approximation for various 
string tension a. The first set is the prediction for only the pion loops, while the second set is with both pion and kaon loops 
included. 



U I LJC V I 










pion loops 






n no 
u.uy 


n 1 m 


-U.loU 


U.U 


n 1 o 
U.lz 


u.uyz 


-U.i4U 


U.U 


0.15 


0.085 


-0.126 


0.0 




pion and kaon loops 






0.09 


0.132 


-0.151 


-0.034 


0.12 


0.121 


-0.133 


-0.032 


0.15 


0.112 


-0.120 


-0.031 



In Tabic IV arc shown the calculated anomalous magnetic moments of the u, d and s quarks for A = 675 MeV for 
various choices of a. Clearly, the results depend on the constituent quark masses. These are given in Table I for the 
considered string tensions. 

Using Eq. (57) the K-term in Eq. (69) yields a nucleon magnetic moment correction 

5/^2 = = 3 {tpsymm \l^g{'i-)(^z{'i-)\ i'symm) Aq (72) 

with 

_ Jr^dr{\g\^-\f\^) 

Jr^dr{\g\^ + m ^'''> 

In Table III the predictions for the nucleon are shown including also the one-pion loop contributions (72) and two-body 

currents. Our results obtained for the one- loop corrections are smaller than reported by Glozman and Riska [39]. This 
is due to the inclusion of the lower component in the single quark orbitals. Neglecting these we recover the results 
of Ref . [39] . Prom Table III we see that the proton and neutron magnetic moment is in reasonable agreement with 
experiment for a string tension of (t = 0.1 GeV"^. For this value of the string tension the model predicts a nucleon 
mass of 940MeV, remarkably close to the empirical value. The anomalous magnetic moment contributions are found 
to be of the order of 10%. 

Due to the one-loop contributions the magnetic moments of the other baryons are modified. Corrections from kaon 
loops have also been considered. Because of the larger kaon mass the contributions are expected in general to be 
smaller in size than those of the pion loops. In Table IV the calculated anomalous moment of the strange quark due 
to the kaon one-loop corrections are given. In the calculations a cutoff mass of A = 675 MeV has been used. The 
isoscalar and isovector anomalous magnetic moment pieces are also changed by the kaon loop contributions. Prom 
Table IV wc sec that the kaon loop contributions are indeed smaller in magnitude as compared to the pion loop ones. 
The full results for the magnetic moments of the baryon octet and decuplet, including the pionic exchange currents 
and the pion and kaon one- loop contributions are summarized in Table V. Por the value of the string tension ct = 0.1 
the overall agreement with the experimental data is reasonable. From the table we see that the anomalous magnetic 
moment contribution leads to an improvement of the predictions. 
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TABLE V. The magnetic moment hb of the haryon octet and decuplet in units of nuclear magneton, including the anomalous 
contribution S^b arising from the pion and kaon one-loop diagrams and the pion exchange corrections for different string tension 
a. Also are shown the experimental results 
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VII. CONCLUSION 



We have written down the general effective quark Lagrangian as obtained from the standard QCD Lagrangian by 
integrating out the gluonic degrees of freedom. Considering the baryon Green's function, neglecting gluon and meson 
exchanges, we find in lowest order of the approximation scheme that it is given by a product of 3 independent single 
quark Green's functions. As a result the Hamiltonian can be written as a sum of three quark terms, where the single 
quark solutions satisfy the Dyson-Schwinger equation with a nonlocal kernel. 

The nonlinear equation for the single quark propagator S (attached to the string in a gauge-invariant way) has been 
solved in the Gaussian correlator approximation. The resulting 3-quark wavefunction has been used to determine the 
magnetic moments of the baryons. This has been done for both the octet and decuplct of the SU(3) flavour group. 

Comparing the predictions we find that the magnetic moments are mostly in close overall agreement with the 
experiment for a string tension of it = 0.1 GeV"^. We find, that the predicted magnetic moment of the nucleon is 
improved substantially once we choose a string tension to give a reasonable nucleon mass. The same applies for the 
A-isobar. Effects due to the presence of virtual mesons are in general expected to be important. We have estimated 
the pionic one-loop and one pion exchange contributions to the magnetic moment. The single quark corrections from 
pionic loops arc foiuid to be of the order of 10%, whereas the total effect of 2-body current contributions are predicted 
to be small, to be contrasted to the results of Ref. [32]. This is due to the cancellation of the pion in flight and 
pair term in the present model. Because of the anomalous magnetic contributions there seems to be somewhat an 
improvement of the predictions. 

Our results for predictions of the magnetic moments of baryons are encouraging, but are in need of including higher 
order corrections. In particular, the mass spectrum obtained from our lowest order approximation docs not contain 
the — A mass splitting. This is due to neglecting contributions like the hyperfine interaction arising from the one 
gluon interaction. It is clearly of interest to investigate how the magnetic moments are changed when effects from 
color Coulomb and hyperfine interaction are accounted for. 
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APPENDIX A: MAGNETIC MOMENT CALCULATION IN COORDINATE SPACE 

The one-quark contribution to the magnetic moment can be written as in (45) 

{AHi) = -e« I (<^(i) V(i) Ax^^) + x^'^ A(^(i))dV, (Al) 

where A = i (H x r) is the vector potential of external constant magnetic field- 
Inserting in (A.l) if^^^ = g{r)fljiM and x^^^ = if{r)^ji'M, and taking into account that flji'M = —{cn)^jiM, one 
easily obtains 

(AiJi) = -ieW J d\{g*f + /*5)r^*jM{('^n)(nH) - t7H}i7,w (A2) 

Eq.(A.2) contains the matrix element / dnfljij^niUki^jiM, which simplifies when / = 0, so that {niUk) = |^ife. 
In this case one obtains, taking into account relation {AHi) = AE = — /x^^^H, 

M^'^ = -l^'^ J i9*{r)f{r) + r{r)9{r))rd'rn*iM<y^'^njiM = 

= _o-(i)^e« j Re{g*{r)f{r)ydr (A3) 

In the case of a local scalar potential U{r) one can further express /(r) through g{r) using the Dirac equation for the 
one-quark state 

rfir) = ^ , i ^igr) + -gr] (A4) 

Introducing (A. 4) into (A. 3) and integrating by parts one obtains 

For U{r) = ar one obtains Eq.(49). 
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APPENDIX B: MAGNETIC MOMENT OF THE MULTIPLET 



In this Appendix the calculation of the nuclcon magnetic moment is generalized to the baryon octet and decuplet. 
By analogy with the fully symmetrical 3g wave function for the nuclcon, Eqs. (46-47), wave functions for the baryon 
multiplets can be formulated. The flavor octet with total spin 1/2 up becomes. 



symm 



symm 



V2 



6 



u+u^d^ — u^u^d^ — u^d^u^ — d^u^u^ + 2'u+'(i-|_(i_} 

— {2u-d+d+ — d^u+d+ ~ u^d-d+ + 2d+u^d^ 
) 

d^d-U+ — d-d+u+ — d+u+d- — u+d+d- + 2d+d+U-} , 



V2 



symm 



{2s-U+u+ — U-S+U+ — S+U-U+ + 2u+S-U+ 



symm 



symm 



— U-U+S+ — u+s+U- — s+u+U- + 2u+u+S-} , 

— T" {u+d_s+ + d+U-S+ + s+d-U+ + s+U-d+ — 2u+s-d+ — 2d+S-U+ 
6 

+U-d+s+ + d-U+s+ — 2s_(i+vi+ — 2.s_7i+fi+ + u-S+d+ + d-s+u+ 
—2u+d+S- — 2d+u+S- + s+d+u- + s+u+d- + u+s+d- + , 
V2 

-— {2s_d+rf+ — d-S+d+ — s+d-d+ + 2d+S-d+ 
6 

—d+d-S+ — d-d+s+ — d+s+d- — s+d+d- + 2d+d+S-} , 

—— {M_rf+s+ — d_u+s+ + U-S+d+ — d_s+u+ — ■u+rf_s+ + rf+u_s+ 



symm 



s+rf_u+ + s+w_d+ + s+d+w_ — s+u+d- — w+s+rf_ + rf+s+u_} , 

+ — u_|_s_s+ + 2s+u_s+ 



V2 , 

— {2u-s+s 
b 



—S+S-U+ — s_s+u+ — s+u+S- — u+s+s_ + 2s+s+u_} , 
\/2 

-— {2d_s_|_s+ — s_d+s+ — d+s_s+ + 2s+(i_s+ 
6 

— s+s_rf+ — s_s+rf+ — s+rf+s_ — rf+s+s_ + 2s+s+(i_} , 
where the subscripts (±) refer to the spin projection. For the flavor decuplet with total spin 3/2 up we have 



symm 



A++ 
symm 

A+ 

symm 

A° 
' symm 



symm 



■ symm 



symm 



symm 



symm 
symm 



1 

v3 

1 

1 

1 
1 
1 



{u+d+s+ + d+u+s+ + u+s+d+ + s+u+d+ + d+s+u+ + s+d+u+} . 
{d+d+s+ + d+s+d+ + s+d+d+} , 

{s+s+c?+ + s+d+s+ + rf+s+s+} , 

S+S+S+. 



(Bl 

(B2 

(B3 
(B4 

(B5 
(B6 

(B7 
(B8 

(B9 
(BIO 

(Bll 

(B12 
(B13 

(B14 

(B15 

(B16 

(B17 

(BIS 



18 



These fully symmetrical wave functions Eqs. (B1)-(B18) can be written symbolically as, 

^V.,™™ = Tt}{hf2f^)i^{^^f^p!;. (B19) 

As the orbital of the s-quark is heavier than the u- and d-quark orbitals Eq. (51) has to be split up in contributions 
from the u, d-quark and from the s-quark. Using the symmetrical wavefunction Eq. (B19) this is realized by writing, 

= 3Mi^) = 3 (rj^7(/i/2/3)V'^ V'^^V'^^ |e,(l)a,(l)| r^^^(A/2/3)^^ V^^V^^) (B20) 
/1/2/3 

with, 

^n=-lj9l{r)ff,{rydr. (B21) 

The flavor index /j can take the values u,d or s. Note that A„ = as the same orbital is taken for the u- and 
d-quark. Evaluating Eq. (B20) for the different baryon wave functions Eqs. (B1-B18) results in the expressions in 
Table VI. 



TABLE VI. The matrix elements of the e.m. current for the haryons. 



Ix 

3 A 

9 At, 

8 \ -1- 1 X 
07 '^u "T 27^3 

2. \ _|_ J_ \ 

— 4 X -I- 1 X 

27^14 -t- 27^3 

^ n As 

-Ax - Ax 

27 As 27'^" 

- 4 \ I 1 \ 
27 As -I- 27 A 



i\ 
1 Au 

1 X 

gAu 



1 \ 

~ 3 Au 
4 X _ 1 X 

9 9 

— 2 \ 1 X 

^Au — ^As 

— Q As + qAi^ 

-2x - lA 
9 ^ 9 ^ 

~3A« 



N 



P 
n 

E+ 
S° 

A 
=•0 



A+ 
AO 
A- 
E+* 

E"* 
E-* 
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APPENDIX C: PIONIC TWO-BODY CONTRIBUTION TO THE MAGNETIC MOMENT 



In this Appendix the pion-in-flight and pair contributions to the magnetic moment of the nucleons are given. 
Following Ref. [33] these contributions are determined by taking the curl of the pionic two-body currents Eq. (68). 
The 3-quark state ^ is given by the product of three single quark orbitals Eq. (35). Because of symmetry considerations 
it suffices to calculate the magnetic moment contribution of pion exchange between say the second and the third quark 
only and multiply the result by a factor of 3 to include the contribution of the other possible permutations of quark 
pairs. 

Considering the pion-in-flight contribution first (Eq. (64)), taking the curl gives rather long expressions which can 
be divided into two parts, 

The first part gives the larger contribution and can be written as, 

[9{k'2)f{k2)m)hkz)^^ (pA - p,A,) 

+9{k'2)f{k2)f{K)9{k3)^^^ ((2p - p' + V3q) A - {2p, - p'^ + V35,)A,) 
+/(fc2)ff(fe)5(A;^)/(fc3)^ ((2p - p' - \/3q)A - {2p,-p'^ - %/3g,)A,) 
+f{k'2rg{k2)mr9{ks):^ (p'A -p'.A,)} x 



1 + 2- 



A2+A2 ; VA"+A2 
The second part comes from the curl applied to the wavefunctions 



[^fik2)mfiks) (k2 X ka)^ - ^fik2)mmks) (k2 X k;)^ 
-^^9{k2)m)fiks) (k; X ks)^ + ^-^9{k2)~m)g{ks) (k; x t;)^ 



+/(fc^)5(fc2)ff(fc^)/(fc3)^ (k2 X ka)^ -~f{k'2)~g{k2)hk'^y9{k^)Y (^2 X k^ 



^2 

,A2 + m2\/ Kl 



A^ + AlJ VA2 + A2 

The normalization factor N is the same as used before in the single quark current contribution (Eq. (61)). The 
momenta are expressed in terms of the Jacobi coordinates Eqs. (62) again, but from imposing the Breit system and 
momentum conservation we now get P' = — P = Q/2 and 2\/3(q' ^ q) = Q- In writing down these expressions use 
has been made of the spin-isospin operator sandwiched between the fully symmetric wavefunctions in spin-isospin and 
orbital space of the 3 quarks, 



(^P (r(i) X t(2)) a^^V^") \=-^V'" (rW X t(2)) a^^V^^) ^A" \ 

\Tsymm V ' ^ 3 r symm / \Tsymm V ' J z i- 3 



symm j (f^^ 



It should be noted that the spin-isospin factor (C4) is identical to that found for the tri-nucleon case. For all the other 
baryon wavefunctions given in Appendix B the matrix element of the considered two-body e.m. operators vanish, 
because of the isospin structure of the e.m. operator. Hence the considered two-body currents contribute only to the 
magnetic moment of the proton and neutron. 
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The second part Sfif is a relativistic effect which enlarges the values by about 10% and which vanishes in the static 
limit as is shown at the end of this section. 

In the same way the pair term can be analyzed. We find 



with 



' ^ 9{k2)f{k2)9{k^)g{kz) 



k2 \3 A2 + m2 

1 /2 , Ak'2-A,(fc^), 



1 /I , Ak3-A,(fc3;. 



1 f2 Ak^ - A,(fc^),^ m)9{k2)f{k',rg{ks)] f t^TT ) (C6) 



and 



/ ^^~g{k2) (3(fc^)/(fc3) (k;k3 - {k'2Uk).) - Rk'srgih) (k2k3 - ik'2)M)^ 



+g{k'srg{ks) [-^r^/(fc2) (k;k2 - (fc^).(fc2).) + ^^^(^2) (k^k^ - ikUkh) 

A? 



- f{k'2rgik2)gik'Mk3)lr {^'2^ - {K)^)^)^ 



,2 X 2 



In the non-relativistic limit the lower component of the wavefunction can be expressed in the upper component as 

/» = -^Kfc), (C8) 

where nig is the constituent mass of the quark. As a result, the pionic two-body current contributions Eqs. (C1-C3) 
and Eqs. (C5-C7) can be simplified considerably. We obtain for the pion-in-flight contribution 



-^M^ = 6r»gfc^ /^'g^W + ~gik'2)g{k2)mr9{k3) 



For the pair term we get 



(C9) 



(CIO) 



while and 5fj,^ vanish. These expressions agree with the results of Refs. [33] and [35]. 
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APPENDIX D: ANOMALOUS MAGNETIC MOMENT CONTRIBUTIONS FROM PION LOOPS 



Our starting point is the e.m. currents, corresponding to the one-loop diagrams shown in Fig. 3, 



/(«) 



(27r) [{p-k) 
Al 



2 — + ie\ [fc^ 
P - ml 



mt 



ie] {k + Q) -ml + ie 



{k + QY - ml 



fc2 -Al{k + Qf _ A2 \ (A: + Qf - A2 fc2 - A2 



and 



/- d*k 

J 



75 (i*'- + 7m (i^- ^ + '^q) 75 



(27r)^ ■ [(p' - fc)2 - to2 + ie] [(p - /e)2 - to2 + ie] [A;2 - m2 + ie] V - A2 



A2 



(Dl) 



(D2) 



Since we have assumed a finite form factor at the nqq vertex, similar as in the two-body current case, the two additional 
terms are needed in the last factor of Eq. (Dl) to satisfy current conservation. ^From these currents the anomalous 
magnetic moment has to be extracted. By applying the Gordon decomposition to the current Eq. (69) near Q'^ = Q 
it can be seen that the anomalous magnetic moment k is the term proportional to —-^Kfj, with isT^ = p'^. To 
isolate this term the currents are rewritten by explicit evaluation of the 7-matrix algebra and taking the limit 0. 
Using the approximation that the initial and final quark is on-mass shell we obtain 



(27r)^ [A;2 - 2pk + ie] [fc2 - ml + ief W - Al 



Al 



1 + 2 



fc2- A2 



(D3) 



and 



■ 2 J 



d^k 



(27r)^ [/c2 - 2pk + ief [fc2 - ml + ie] V - 



A 



As the tensors C^" depend only on the initial and final momenta they can be written as, 



(D4) 



(D5) 



where An'' are Lorentz invariants. It can readily be seen, that only the first term A^*-* contributes to the magnetic 
moment. Substituting Eq. (D5) in Eqs. (D3-D4) and taking the initial and final quark on-mass shell we find for the 
anomalous magnetic moment corrections 



^1 I 



(D6) 



-4iM 



pmqg: 



irqq 



1 A.ib) 
^ ^1 • 



The Lorentz invariant expression vl^*^ can immediately be determined from the tensor C^^ . We get 



(i) 



At^ = 



(D7) 



(D8) 



Inserting Eq. (D8) in (D6-D7) the expressions (70-71) are obtained. 

The kaon one-loop diagrams can be calculated in similar way. The starting point is the expressions Eqs. (Dl- 
D2) again, where the mass of the pion is replaced by the mass of the kaon and the isospin structure is changed to 
{tz -|-3Y')/2 and — (| -|- |F) respectively in Eqs. (D1-D2) with Y the hypercharge. The expressions for the anomalous 
magnetic moment due to the kaon loop become, 
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(D9) 



and 



(6)_-2 /2 4 \ 2Afp r d^fc 4(p-fc)^-p^fc^ + 3m,(M,-m,)p-fc / A 



Q 



y2 



with Mg the mass of the intermediate quark, the mass of the initial and final quark. The coupling constant QKqq 
and the cutoff A^ are taken the same as for the pion loop. 
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